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Abstract. 

Gravitational lensing has developed into one of the most powerful tools for the analysis of the 
dark universe. This review summarises the theory of gravitational lensing, its main current 
applications and representative results achieved so far. It has two parts. In the first, starting 
from the equation of geodesic deviation, the equations of thin and extended gravitational 
lensing are derived. In the second, gravitational lensing by stars and planets, galaxies, galaxy 
clusters and large-scale structures is discussed and summarised. 



1. From general relativity to gravitational lensing 

1.1. Preliminaries 

We assume throughout that a valid model for the geometry and the evolution of the universe 
is given by a Friedmann-Lemaitre-Robertson-Walker model with cosmological parameters 
narrowly constrained by numerous cosmological observations. The observational evidence 
for this cosmological standard model has recently been reviewed elsewhere [fl6l . Here, we 
just summarise its main parameters in Tab.[T] adapted from H262L 



Table 1. Main cosmological parameters, adapted from [262 ] . The Hubble constant h is 
dimension-less and defined by Ho = lOO/zkms -1 Mpc -1 . "CMB only" means parameters 
obtained from the analysis of the 7-year WMAP cosmic microwave background data, "CMB, 
BAO and Hq" takes additional constraints from baryonic acoustic oscillations and external 
measurements of the Hubble constant into account. See 12621 for detail. 



Parameter 


CMB only 


CMB, BAO and H 


Hubble constant 


h 


0.710 


+ 


0.025 


0.704 


+ 0.013 
0.014 


baryon density parameter 




0.0449 


+ 


0.0028 


0.0456 


+ 0.0016 


cold dark matter density 


^cdm.O 


0.222 


+ 


0.026 


0.227 


± 0.014 


cosmological constant 




0.734 


+ 


0.029 


0.728 


+ 0.015 
0.016 


power-spectrum normalisation 


0"8 


0.801 


+ 


0.030 


0.809 


+ 0.024 
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1.2. Equation of geodesic deviation 

Gravitational lensing studies the effects of light deflection on the appearance of cosmic 
objects. Many approaches to gravitational lensing exist, some of them heuristic, others rather 
formal mathematically; see Il423l 14921 13581 14661 l425ll for examples. 

We discuss gravitational lensing here under three main assumptions which are underlying 
the entire treatment. First, we shall remain within the framework of general relativity. Second, 
we shall assume that gravitationally lensing matter inhomogeneities have weak gravitational 
fields in the sense that their Newtonian gravitational potential is small, O «; c 2 . Third, the 
sources of the potential are assumed to move slowly with respect to the mean cosmic flow, 
such that peculiar velocities are small compared to the speed of light. The assumption of 
weak, slowly moving gravitational lenses is well valid in all astrophysical applications except 
for light propagation near compact objects, which is not covered by this review. Lensing by 
moving and rotating astrophysical bodies has been discussed in the literature and generally 
been found to be negligibly small flM I23U 14431 [T87l I4T7H . Within the frame of these 
assumptions, gravitational lensing can be considered as a complete theory whose theoretical 
aspects are fully developed, including the mathematics of singularities in the lens mapping 
(cf. the textbook by IT388H ). 

Our approach here begins with a congruence of null geodesies (see Fig.[T]), modelling a 
bundle of light rays in the approximation of geometrical optics within general relativity. We 
select one of the light rays as a fiducial ray and parameterise it by an affine parameter A which 
we shall later specify (cf. Fig. [T]). The tangent vector to the fiducial ray is 
dx^ 

k M = , (1) 
dA 

and we choose to normalise k such that its projection on the four-velocity u obs of a freely 
falling observer is unity, 

(k, u obs ) = 1 . (2) 

Since the wave vector k of a light ray, when projected on w obs , gives the (negative^]) frequency 
o» obs measured by the observer, we must choose A such that 

~ k = 7^ = W L ^r (3) 

When the fiducial ray is connected with a neighbouring ray by a curve y(<x) with tangent 
vector da-y = v, the equation of geodesic deviation or Jacobi equation 

V\v = R{k, v)k (4) 

quantifies how the vector v changes along the fiducial ray in response to the curvature. 

To evaluate this equation, we introduce a two-dimensional screen perpendicular to the 
four-velocity u obs of the freely-falling observer and to the normalised tangent vector k of the 
fiducial light ray, i.e. a screen in the 3-space of the observer perpendicular to the light ray. 

X We adopt the signature -, +, +, + for the metric. 
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Figure 1. A congruence of null geodesies with the fiducial light ray, a curve y connecting rays 
in the bundle, and the tangent vectors to y and the fiducial ray. 



This screen is spanned by the vectors E\£, which are parallel-transported along the fiducial 
ray, 

Vjf^i = . (5) 
Denoting with v 1 ' 2 the components of v on this screen, the equation of geodesic deviation turns 



into 



Vi. 

k\ „2 



= T 



H436B . where the optical tidal matrix 

3(T) n- %(T) 



T = 

with the components 



(6) 



(7) 



1 



1 



H = — -R a pk c K + -C a p y se a 



y *6 



(8) 



appears. It contains the Weyl curvature specified in Eq. (35) below, and the complex vector 
6 = E\ + iE 2 . 
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We continue with the fundamental assumption that we can split the optical tidal matrix 
into two contributions, 

T = T bg + T d , (9) 

where the first is due to the background, homogeneous and isotropic space-time 
of a Friedmann-Lemaitre-Robertson- Walker model, and the second is due to local 
inhomogeneities described by their Newtonian gravitational potential. This assumption is 
justified by our fundamental assumptions laid out in the beginning, i.e. that the peculiar 
velocities of any inhomogeneities with respect to the comoving flow of the background matter 
are small compared to the light speed, and that their gravitational fields are weak in the sense 
that their Newtonian potential O is small, O/c 2 =: tp «. 1. 



1.3. Background contribution 

We first work out the background contribution T bg to the tidal matrix. According to 
our assumption, the background has the Friedmann-Lemaitre-Robertson- Walker metric with 
spatial curvature K, 

ds 2 = a 2 (r]) [-dr] 2 + dw 2 + fl(w)dQ. 2 ] , (10) 

where r\ is the conformal time related to the cosmic time t by adr\ = cdt, and /k(w) is the 
comoving angular-diameter distance as a function of the comoving radial distance w, 

I -4= sin ( V^w) K > 
fn(w) = \ w K = . (11) 

I -±= smh(^^Kw) K < 

In the homogeneous contribution to the optical tidal matrix, the Weyl curvature must vanish 
because of the symmetry of the background model. Since k is a null vector, we can further 
write 

R a ik a F = G aj3 k a W = ^-T a/3 ~k a F . (12) 

Here, Einstein's field equations enter for the first time. In alternative, metric theories 
of gravity, the relation between the global properties of light propagation and the matter- 
energy content of the universe will change, but the description of light deflection by local 
inhomogeneities in weak gravitational fields remains valid to some degree, as specified in 
Subsect. |1.5| below. 

Inserting the energy-momentum tensor of an ideal fluid and assuming negligible 
pressure, p = 0, we have 

T ap k a ¥ = pc 2 (u,k) 2 . (13) 

Now, the projection of k on u is the frequency of the light as measured by an observer co- 
moving with the fluid. Since we normalise k according to Eq. ([3]) by the frequency a> bs 
measured by us as observers, we must have 

\(k, u)\ co 

\(k,u)\ = 11^1 = = i +z (14) 

^obs ^obs 
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according to Eq. ([3]), where z is the redshift of the fluid with respect to the observer. For a 
pressure-less fluid, p = p (l + z) 3 with the density p measured by the observer, and thus 

K = -^ r p Q (l+z) 5 . (15) 

c 1 

The optical tidal matrix of the homogeneous and isotropic background is thus 

r bg = -ni 2 , (16) 

where J 2 is the unit matrix in two dimensions. 



To satisfy Eqs. (|1| and ([14]), we now have to choose A such that \(k, u)\ = 1 +z. Assuming 
that peculiar velocities can be neglected and that we can thus project k on the mean flow 
velocity = 6^, we must have 

,r v l dx \ dx ° cdt 1 -1 

*"> = («•")= dT = dI = 1+z = a • (17) 

and thus dA = acdt = a 2 drj. The equation of geodesic deviation for the smooth background 
then reads 

v ^ = d5 = ^> ;=w ' (18) 

where we have used in the first equality that the basis vectors E t spanning the screen are 
parallel-transported along the fiducial ray, see Eq. ([5]). 

We now replace dA by dw, the comoving radial distance element, using that 

dw = drj = a~ 2 dA => dA = a 2 dw (19) 



for light, since ds = in Eq. ( 10). Introducing further the co-moving bundle dimensions v'/a, 
their propagation with w is given by 

i_ ( -) = a 2 — ( v"a - v'a') = a 2 ( v u 'a - v'a") , (20) 
d«r \ a ) dA y ' v ' 

where the primes denote derivatives with respect to A. Moreover, since dA = acdt, we have 

da 1 da a 

^ = tt = --77 = - (21) 
dA ca dt ca 



da' 1 da da' 1 a d a Id /a\ 2 



and 

— - — — ~' ~" - — |~| (22) 

dA ac dt da c 2 a da a 2c 2 da \ a i 
We can now insert Friedmann's equation 

<a\ 2 _SnGp A Kc 2 
\al 3 a 3 3c 2 a 
for (a /a) 2 to find 

4nG po K 

c 2 a 4 a 

Note that the cosmological-constant term drops out. Loosely related in this context, arguments 
that the cosmological constant might measurably affect gravitational lensing by isolated 



a" = r - t + - 3 ■ (24) 
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objects [I407H were refuted by 1145511 . The propagation equation (|20|) for the comoving bundle 



dimensions in the homogeneous and isotropic universe thus becomes 
d2 I A 3 u, . 4nGp t T J T J 



-, I , oV' + — r ^v' -K- = -K- (25) 
dw 1 \a c L a 1 a a 



where Eqs. ( 15 ) and ( 18 ) were inserted in the last step, or 
d 2 \ if 



dw 2 



+ K\-=0. (26) 



a 



1.4. Clump contribution 

In presence of local perturbations characterised by the Newtonian gravitational potential 
O = 4>c 2 «; c 2 , we can write the perturbed Friedmann-Lemaitre-Robertson- Walker line 
element as 

ds 2 = a 2 (r]) [-(1 + 20)d?7 2 + (1 - 20) (dw 2 + f 2 (w)dQ. 2 )] (27) 
if the perturbations move slowly with respect to the mean cosmic flow (cf. the more general 



remarks on light propagation in weakly perturbed metrics in the following Subsect. 1.5 ). Our 
next assumption is that the perturbations are well localised, i.e. that their spatial extent is 
much smaller than the curvature scale of the background universe. Then, we can consider 
the perturbation as being locally embedded into flat space and approximate /#(«;) « w. 
Consequently, we shall first study light propagation in the comoving Newtonian metric 

ds 2 = -(1 + 20)d7/ 2 + (1 - 2cf>)dw 2 (28) 

and later combine it with light propagation in the background universe. 
The line element Eq. ( [28] ) suggests introducing the dual basis 

8° = (1 + 0)d?7 , ff = (1 - (p)dui , (29) 

in which the metric becomes Minkowskian. The following calculation is linear in the sense 
that terms of higher than first order in (p are neglected, as are time derivatives of the potential. 
By Cartan's first structure equation, the connection forms are 

w° = , J j = -<Pjff + $0 , (30) 

where </> 7 = dj<p and partial derivatives are taken with respect to comoving Cartesian 
coordinates. Cartan's second structure equation gives the curvature forms 

Q ,. = $ ik & A 6° , = -4> jk & A ff + <f> ik & A 6 j , (31) 

which yield the only non-vanishing elements 

RoiOj = <Pij - Rijij = <Pii + <Pjj : 

^1213 = 023 5 ^1223 = — 013 5 ^1323 = 012 (32) 

of the Riemann tensor. The Ricci tensor and the Ricci scalar are 

R a p = V 2 (f>I A , R = 2V 2 , (33) 
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respectively, where J 4 is the unit matrix in four dimensions. The Einstein tensor is 

G a/3 = V 2 0<^, (34) 

with V = o'di. The Weyl curvature, which we need for the clump contribution to the optical 
tidal matrix, is 

CafiyS = RafrS ~ Ga[yRs]p + Gp[yRs]a + -^9a[yGs\li , (35) 

and has the only non-vanishing components 

1 ->2 2 ->2 

CoiOj = 4>ij - ^ V 077,7 ' C W = 0" + 077 ~ 3 " ' 

Cl213 = 023 5 C1223 = _ 013 > C1323 = 012 (36) 

here. We can now evaluate the clump contribution first in the comoving metric given by 



Eq. ( [28] ) . With the components Eq. ( [32] ) of the Ricci and Eq. ( |36[ ) of the Weyl tensors, we find 
ftd = -V 2 , r d = - (011 - 022) - 2i0i 2 , (37) 
hence the clump contribution to the optical tidal matrix becomes 

T * = ~ 2 ( t U t n ) > = - 2 W • (38) 



012 022 

Let now (O) be the gravitational potential passed by the fiducial ray. Then, the potential 
gradient can be expanded as 

d i (0 - (O) ) = = djd i( f> \ q x j = - 1 - (T d )ij x j . (39) 

We can thus bring the equation of geodesic deviation for the clump contribution into the form 

d 2 x' d 2 y 

^ = ^-™*- m 

where the first equality follows from dA = dw in the local approximation Eq. ( [28] ) to the metric. 
Since the absolute value of the potential is of no importance, we can rename the potential 
difference 6(f> to the potential 0. Finally, we must be aware that the bundle dimensions x l 
evaluated in the local frame are comoving bundle dimensions v'/a in the cosmological frame. 
If we now combine the global, homogeneous and the local, clump contributions 



expressed by Eqs. ( [40] ) and ( [26] ), we find the inhomogeneous propagation equation for the 
comoving bundle dimensions x' 

dw 2 



+ K\x' = -2d l 4> . (41) 



1.5. Remarks on light propagation in other metric theories of gravity 

At this point, it is worth noticing that part of the results derived so far remain valid beyond 
general relativity. In a metric theory of gravity, the most general, lowest-order deviation from 
the Minkowski space-time is 

ds 2 = -(1 + 20)d77 2 + (1 + 20)dw 2 (42) 
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with the two otherwise unspecified Bardeen potentials <p and if/ with 0, if/ «: 1. Since null 
geodesies are conformally invariant, light propagation can be described by the null geodesies 
of the conformally transformed metric 

ds 2 = (1 + 20) [-dr] 2 + (1 + 2(A)dw 2 ] , (43) 



with if/ = if/-(f>. The geodesic equation of the conformally transformed metric (43 ) is identical 



to (40) with -2(p replaced by if/. Locally, light deflection is thus governed by the difference 
between the Bardeen potentials which, in general relativity in absence of anisotropic stresses, 
is if/ = -2(f>. The gravitational field equations only enter through the relation between the 
Bardeen potentials and the matter-energy content of the space-time, which also governs the 
global, cosmological properties of light propagation. Thus, in any metric theory of gravity 
which locally contains special relativity, light deflection by localised, weak gravitational fields 
is determined by if/. Alternative field equations will change global light propagation and the 
way how matter and energy source if/. 

1.6. The lens mapping and the lensing potential 
The Greens function G(w, «/) of the operator d 2 /dw 2 + K is 
G(w, w) = sin [ VK(w - «/)] @(w - w) 

= f K (w - w')@(w - uf) , (44) 
where we have introduced the comoving angular diameter distance /k(w) of the Friedmann 



metric, Eq. ( 10). The appropriate boundary conditions for the solution of Eq. (41 ) are 



dx l 

= 0, — =ff, (45) 
w=0 dw w =o 

because both rays start at the observer and enclose an angle (0 l ,6 2 ) there. Thus, the solution 
is 

x\w) = f K {w)ff - 2 I dw'f K (w - w')ff(f> (x j (w'\ w) . (46) 

<J 

In line with the Newtonian approximation, (p c 2 , we assume small deflection angles 
and use Born's approximation to integrate along the unperturbed light path x'(w') x f K (w')8 l in 
Eq. ( [46} . Since numerical simulations show that Born's approximation is typically very well 
satisfied [11061 14441 11881 . we shall adopt it throughout. See [4TJ for a formidable calculation 
to second order. Integrating from the observer to a source at the comoving radial distance w s 
and defining x'(w s ) = fK(w s )f3', we find 

fit = ff-2 P dw' (f K (w')0j, uf) . (47) 
Jo Ik(w s ) v ' 

We define the reduced deflection angle by 

«W):=2 P dw' fK( \) - "° ff<p (friw'W, uf) (48) 

Jo flcM V I 
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and obtain the lens equation 

= ff - a \e j ) . (49) 
Introducing derivatives with respect to angular coordinates 9 j on the sky by 

d x = ~jT7~\ . (50) 
Jk(w) 

we can write the reduced deflection angle as the angular gradient of an effective lensing 
potential 

He j ) = 2 r s dw' ^ w ;-f\ <p(f K (w')e\w') . (51) 

Jo Mu)')f K (w s ) v ' 
The derivatives with respect to 6 have to be interpreted as covariant derivatives on the sphere 
if the curvature of the sphere becomes important. Then, different suitable and convenient 
bases are used. We shall later return to this issue. For the purposes of gravitational lensing, 
the flat-sky approximation is often, but not always justified. Unless stated otherwise, partial 
derivatives are taken with respect to angular coordinates from now on. 



1.7. Local properties of the lens mapping 



Having arrived at this point, it may be useful for clarification to describe lensing without 
reference to specific coordinates on the sphere. The lensing potential if/ assigns a number to 
each point on the observer's sky, 

ifr : S 2 -» R , p i-» if/(p) . (52) 

The gradient of \p defines a vector field a on S 2 by 

a b = dtff , (53) 

where a b is the 1-form dual to a. The codifferential of a b is the Laplace-de-Rham operator of 
if/ and defined to be twice the so-called convergence 

6a b = (*d*d) = 2k . (54) 

The lens mapping defines a map (p from the observer's sky to a sphere on which the sources 
are located, 

<p : S 2 -» S 2 , p h-> (f(p) . (55) 

Its differential Dip describes how sources are locally deformed under the lens mapping. In 
absence of scattering, absorption or emission, the phase-space distribution function / satisfies 
Liouville's theorem. This implies that / oc a)~ 3 I(a>) is constant along null geodesies, where to 
and I(co) are the frequency and the specific intensity of the light. If the frequency is unchanged 
by the lensing mass distribution, I{oS) is constant and the flux from the source is changed only 
because the lens mapping changes the solid angle under which the source appears. Thus, 
lensing causes the magnification 

H = IdetfZty)!" 1 ■ (56) 
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Points p cr it where the lens mapping ip is singular, det(Dip)(p CT it) = 0, are called critical points. 
They form closed curves, the so-called critical curves. Their images <p(p C nt) under the lens 
mapping are called caustics. 

If we now introduce charts hi : S 2 D U — > C/' c K 2 with p e U and /i 2 : S 2 D V — » 
V c K 2 with <£>(p) 6 V, we can locally span U' and V by Cartesian coordinates (#',# 2 ) and 
(J3 l ,j3p), respectively, and express the Jacobian matrix Dtp by 

OB 1 da 1 ■ d 2 d/ 

W'-£-4-S7-4-s£»- :4 >-*i- 

It is convenient and instructive to separate Dip into its trace and a trace-free part, 



where 



= ~<9 ! <9^ (59) 



is the convergence introduced as the codifferential (divergence) of a in Eq. (54), and 

71 = 2(^1-^2). 72 = ^2 = ^1 (60) 
are the components of the shear 

7 = 7i+i72- (61) 

While the convergence is responsible for stretching a source isotropically under the lens 
mapping, the shear is responsible for its distortion (see Fig. [2]). In fact, a circular source 
of unit radius is mapped into an elliptical image with semi-major and semi-minor axes 

A = (\- K -\y\r l , B = (l-^ + | 7 |)- 1 . (62) 

This distortion allows the systematic detection of gravitational lensing. The relative axis ratio 
of elliptically distorted images, the ellipticity e, is the so-called reduced shear 

e = X+b = T^k = 9 ' (63) 

This illustrates that not the shear y can be directly measured from image distortions, but only 
the reduced shear g. This is because the absolute size of a source is typically unknown and 
image ellipticities are invariant under transformations 

Dip -> A(Dip) with A ± , (64) 

because they change only the size of the image but not its shape and leave the reduced shear 
unchanged. The convergence is transformed by Eq. (64) as 

k-*1-A(1-k). (65) 

With A = 1 - 6k and 8k <sc 1, this transformation is approximately 

k — > k + 8k (66) 

if k «; 1. Then, it corresponds to adding a sheet of constant surface-mass density 8k to the 
lens, whence it has been called the mass- sheet degeneracy II 1201 11531 . 
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Figure 2. The effect of shear y (left column), first flexion F (middle column) and second 
flexion G (right column) on an idealised circular image. 



In the coordinates introduced by the charts h\ 2 before Eq. (57 ), the magnification is 

(i = \det(D<pT l = ^—^ - . (67) 



(\-Kf-y\-y\ ' 
In these coordinates, the critical curves are given by 

{\-K) 2 -y\-y 2 2 =Q. 



(68) 



Third-order derivatives of the lensing potential, termed gravitational flexion [115 1 L |5), 
may become measurable in the near future. Since partial derivatives of commute, only four 
of the original eight components of djdjdtiff = ifrijk are independent, e.g. 0m, 0m> 0m and 
if/222- They are commonly combined into the first flexion F with the two components 



^i = 2 ('Am + ^122) » 



and the second flexion G with 
1 



G\ = - (^111 



30i 22 ) , G 2 = - (30112 - 0222) ■ 



(69) 



(70) 



Applying the local definition of the convergence in Eq. (59) to the lensing potential 
Eq. (151}, using Eq. (|50j>, we find 



k[x ] ) = dw — — d x ,d <p[x ] ,w , 

v ; Jo Mw s ) v > 



(71) 



where the partial derivatives of <p under the integral are to be taken with respect to the co- 
moving coordinates x l perpendicular to the line-of-sight. We can, however, augment this 
two-dimensional Laplacian by its component along the line-of-sight, 



d x ,dU = v 2 = 4 -^ P 



(72) 
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because the derivative with respect to the line-of-sight direction vanishes after the line-of- 



sight integration in Eq. ( |7Tj ) if the potential is localised. Inserting Eq. ( [72] ) into Eq. ( |7Tj ) shows 
that the convergence is the suitably scaled and geometrically weighed surface-mass density of 
the lensing matter inhomogeneities, 

K(6) = — r dw — — p \f K (w )ff, w ) . (73) 

c 2 Jo MH) v ; 

2. Thin gravitational lenses 

2.7. Introduction and axially -symmetric lenses 

In many astrophysically relevant applications of gravitational lensing, the lensing mass 
distribution is very thin compared to the cosmological distances between observer, lens, and 
source. We can then introduce the thin-lens approximation, replacing the Newtonian potential 
4> by its line-of-sight projection 



(p -> S(w - w d ) J dw'<f> =:6(w- w d )((> (Z> , (74) 

located at the co-moving radial distance w d from the observer. From now on, unless stated 
otherwise, we shall assume that the spatial curvature vanishes, k = 0, allowing us to identify 



A<2) 



fk(w) = w. Then, the lens mapping Eq. (47) becomes 



; = Qi _ 2 W d (w s -W d ) Qi (2) 

w s v ; 

where the partial derivative is again to be taken with respect to the angular coordinates ff ' . The 



lensing potential Eq. (51 1 simplifies to 



^ff) = 2^^^ (w d ff) , (76) 



and the deflection angle simplifies accordingly. The convergence Eq. (73) becomes 
accordingly 

• 4nG (w s - w d )w d C , i • a 
k{6 ) = — dw p [w 6 , w J . (77) 

The line-of-sight projection of the mass density p is the surface-mass density 

2(6/') = J dw'p (w'ff, w) , (78) 

and the quantity 2 cr defined by 

, AnG (w s - w d )w d 

c 1 w s 

is called the critical surface mass density. With these definitions, the convergence of a thin 
lens is 

k(0) = ^ . (80) 

The light-travel time along the deflected light path is longer than along the undeflected 
path for two reasons. First, the deflected path is geometrically longer, and second, there is 
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source 



observer 



w. 



ds 



lens, 
deflector 



angle difference 0-/3 



Figure 3. Illustration of the co-moving undeflected and deflected light paths and their length 



difference (see Eq. ( 82 1 



a gravitational time delay due to the local gravitational potential causing the deflection. The 
geometrical time delay is most easily worked out in the conformally static, spatially flat metric 

ds 2 = -drf + dw 2 , (81) 

but the following derivation can be straightforwardly generalised to the spatially curved cases. 
We construct a triangle between the source, the deflection point and the observer and introduce 
its side lengths w d from the observer to the deflection point, w ds from the deflection point to 
the source, and w s from the observer to the source (see Fig. [3]). From the cosine law, we know 
that 

w ds = if s + w d - 2w s w d cos(9 — P) ~ (w s - u>d) + w s w d (9 -J3) , (82) 

where the small-angle approximation for the cosine was used in the last step. In the same 
approximation, we can rewrite this equation as 

w ds *ws-w d + ™ sWd (0 - pf (83) 
2(w s - w d ) 



and invert it to approximate 



w s - w d ~ w ds - w d w s (9 - P) 2 (84) 

at the same order in (9~P). The co-moving geometrical difference between the deflected and 
the undeflected path lengths is thus 

Aw = w d + w ds - w s » — -{9 -p) « . (85) 

2(w s - Wd) w ds 2 
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Since the conformal time delay A77 = Aw; is very small compared to the Hubble time, we can 
write the cosmological time delay Af geo = Arj/c or 



Af Reo (0) = — 



Aw w s w d (6 - f3) 2 



c cw d& 2 

At the same time, the light experiences the gravitational time delay [144511 

(f>(6w',w') 



Afgrv(0) 



= -2 J dw' 



(86) 



(87) 



which, when expressed in terms of the lensing potential of a thin lens (76 ), can be written as 

1 W s lf d 
C W d; 



At g „(e) = 



(88) 



The combination of both the geometrical and the gravitational time delays gives the total time 
delay 



At{9) = - 

c w d 



1 w s w d ((9-J3) 2 



(89) 



The stationary points of At, where d t At = 0, satisfy the lens equation. This illustrates Fermat's 
principle in gravitational lensing: Images are formed where the light-travel time is stationary. 
As an immediate corollary, this implies the odd-number theorem: non-singular thin lenses 
produce an odd number of images. In absence of the lens, At is a paraboloid in 9 for fixed /3, 
which has a single minimum where one image is formed. Further stationary points are added 
in pairs as local maxima and saddle points. The total number of images must thus be odd. 

A further important result is that the light-travel time scales like an effective distance 
w s w d /w dfi , which is itself proportional to the Hubble distance c/H . The time delay between 
different gravitationally lensed images is thus proportional to the Hubble time Hq 1 . This 
allows the Hubble constant to be derived from measured time delays in multiple-image 
gravitational-lens systems, if ifr and are known, e.g. from lens modelling [140211 . 

Models for thin gravitational lenses are typically constructed by fixing either the effective 
lensing potential iff or the convergence k. Since the two are related by the Poisson equation 
2k = V 1]/ in two dimensions, the lensing potential can easily be obtained by means of the 
Greens function of the Laplacian in two dimensions, 



-II 



d 2 6' k{9') In 



6 -9' 



(90) 



The situation simplifies considerably for axially symmetric lens models. In this case, the 



lensing potential from Eq. (|90|) is 



m = 2 

giving the deflection angle 
a(0) = if/' = 



r*0 s~*oo 

InO I ( 0Wln<9V(#') 

Jo Je 



m{6) 



with m(8) = 2 



f e'de'K{e') 

Jo 



(91) 



(92) 
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Since m is an area integral over the scaled surface mass density of the lens, it is proportional 
to the lens' mass. The primes denote derivatives with respect to the angular radius 6. Using 
the definition of the convergence, 

we can write the components of the shear 

cos 20 / sin 20 / 0'\ 

7i = — =— ^ - ~Z > T2 = — =— 7T ( 94 ) 



2 r 9 ] ' 2 V 6 

in the form 



yi = cos 20 (/c- , 72 = sin 20^- . (95) 



The Jacobi determinant 



det(£V) = (1 - k? - y\ - y\ = (1 - k? - (k - = (l - 2k + fj (l - |) (96) 

can be brought into the form 

/ d m{6)\( m(6)\ 
det(D^) = 1 - --^ 1 - (97) 



de e )\ e 2 

by means of Eq. ([92]). Axially symmetric lenses thus have critical curves where either the 
slope of enclosed dimension-less mass m(8) divided by 6 equals unity, 

or where the enclosed mass equals 6 2 , 

m(0) = 6 2 . (99) 

Because of the preferred directions of distortion occuring there, critical curves described by 
Eqs. (98 ) and ( [99] ) are called radial and tangential, respectively. 



2.2. Simple lens models 

We begin with the instructive case of a point mass M at the coordinate origin, whose mass 
density is 

p(x) = M6(x) . (100) 

Thus, its convergence is 

-> AnGM (We — Wa)wa -> 
K0) = ^ — — 6( Wi ff) , (101) 

and the lensing potential turns out to be 

^ AGM w* — Wa -> 
H0) = — -ln|0|. (102) 

Together with other lensing potentials, this is shown in Fig. [4] With the deflection angle 

, , AGM w s - w d ff 

a = d\p = — — , (103) 

c z WdW s tr 
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the lens equation turns into 

AGM w s -w d 

P = 0-—2Z-- ' (104) 

which is effectively one-dimensional because of the axisymmetry of the lens. Since this is a 
quadratic equation for 6, it always has two solutions unless the source is placed exactly behind 
the lens, i.e. at = 0. Such sources are mapped into rings with radius 

\ C L w d w s 



the so-called Einstein radius. With this definition, the lens equation ( 104 ) can be written 

y = x--, (106) 
x 

where y = fi/6 E and x = 6/6 E . It has the two solutions 

x ± = l -(y± V?/ 2 + 4) . (107) 
The determinant of the Jacobi matrix is 

det(£V) = 1 - -r , (108) 

XT 

thus the critical curve is a circle with radius x = 1, while the caustic degenerates to the point 
y = 0. The two images at x± are magnified by 

^ = -^=4. (109) 

hence their total magnification is given by 

y 2 + 2 

fi = \ M+ \ + \M-\= V- (HO) 

yyjy 2 + 4 

This point-mass model and collections thereof are commonly used for studying the 
gravitational lensing effects of stars and planets. 

A model frequently used for a simple description of galaxy-scale lenses starts from the 
density profile 

Pto = T77-4t2 ' < m ) 
2nG r £ + r£ 

which is characterised by a velocity dispersion <r„ and a core radius r c . This density profile is 
called the non-singular isothermal sphere. It reproduces the approximately fiat rotation curves 
of stars orbiting in galaxies and is thus considered appropriate for modelling gravitational 
lensing by galaxies. Its surface mass density is 

cr 2 1 

1(0) = —^ -, (112) 

2Gw d J0[T(P 

where 9 C = r c /w d is the angular core radius. The convergence can be written in the form 
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The effective lensing potential, 



4 



can be approximated by 

for 9 » # c . The deflection angle 
2*o 



6> 2 + Ql - 6 C In 



e c + ^|^T^\ 



<*(0) 







(114) 



(115) 



(116) 



becomes constant for 6 C — > 0. In the same limit, the singular isothermal sphere has two images 



at 



(117) 



6 ± = ± 2k q 

if/3< 2/f and a single image at 6+ otherwise. 

Numerical simulations show that density profiles of dark-matter haloes are well 
approximated by 



p(r) = 



x := 



x(l + x) 2 ' r s 
This density profile has the convergence 0~5 



2k. 



x 2 -l 



1 - 



^fT^x 2 



arctanh 



1 -x 
1 + x 



and the lensing potential 
HO) = 4* s 



-In 2 
2 2 



2arctanh 



1 -x 
1 + x 



*, : = 



(118) 



(119) 



(120) 



(121) 



(122) 



Peculiar cases are the matter sheet of constant convergence * , whose potential is 

m = K jo 2 , 

and the potential of a constant shear (71,72) with vanishing convergence * = 0, 

«K0) = y(0i-02) + 720102- 

Concluding this brief discussion of axially-symmetric, simple lens models, it should 
be emphasised that these should be seen as building blocks of more realistic and complete 
lens models. Often, real lens systems can be modelled only by embedding simple models 
into external shear, e.g. by adding a potential of the form Eq. ( |122[ ), by elliptical distortions 
of formerly symmetric models, or by combining multiple simple components into one lens 
system. Note that, since we are working in linearised gravity, the superposition principle 
holds and the potentials of individual lenses may be added. 
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Figure 4. Radial dependence of three axially-symmetric lensing potentials introduced here. 



2.3. Reconstruction of mass distributions 

Localised lenses can usually be well described in the flat-sky approximation, i.e. on tangent 
planes to the sky. In such planes, we can expand the lensing potential into a Fourier series and 



write the relations Eq. (59) and Eq. (60) in the algebraic form 



I 2 . I 2 -I 2 . 

k = - 2& ' h = — ' ^ 2 = ■ ( 123 ) 

Hats above symbols denote Fourier transforms, and I is the two-dimensional wave vector. The 
shear components j\ and j2 can be measured from the distortions of lensed images, while the 
scaled surface mass density k is to be recovered. The second and third relations ( 123[ ) yield 
the two independent estimators 

I 2 . . I 2 



* 1 = jp3*' k2 = W 2 h (124) 

for the convergence. In the linear combination 



Kest = Ml + (1 - Wl , (125) 

the weight A can now be chosen such that the variance 

O-l = ((ftest - <*» 2 ) (126) 

is minimised with respect to A. Setting 

0, (127) 



OA 

using that the shear y should be isotropic on average, (y 2 ) = (y\), that its components should 
be independent, {y^) = 0, and that the average convergence should vanish, (k) = 0, gives 
the minimum- variance weights 

(JR - 1 2 2 ) 2 2l 2 l 2 
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and thus the minimum-variance convergence estimator 

ft = y l + ^hlly 2 = C os lifh + sin 2<py 2 = % (e^y) , (129) 

where ip is the polar angle of the wave vector /, and we refer to Eq. ( |o"lj ) for the definition of 
the complex shear. Therefore, in Fourier space, the convergence is obtained from the shear by 
multiplication with a phase factor. In real space, this corresponds to the convolution 

/-> -> -> o -> 2 ~ ft? + 2ift 02 
D(0 - ft)y(ft) d 2 0' , D(0) = — 2 — . (130) 

This opens one important approach to mass reconstruction H237L 



3. Cosmological weak lensing 

3.1. Weak-lensing power spectra 



The lensing potential given in Eq. (51 ) is the fundamental quantity for studying gravitational 
lensing beyond the thin-lens approximaton. We first expand the Newtonian gravitational 
potential into three-dimensional Fourier modes with respect to co-moving spatial coordinates 

x, 

l3. 



<^v)= f ^kq)^*, (131) 



r d 3 < 

J (2* 

introduce the power spectrum 

(kfifrifj) = (2nf8{4 ~ <f)PM) , (132) 
and find the angular correlation function 

fd» r f gLPM*™ (i3 3 ) 

\ ' Jo Jo w,w w & w' J (2ny 



for the effective lensing potential. The average in Eq. (132) is taken over a ll Fourier modes 
with fixed wave number q contained in the accessible volume. In Eq. (133), x = (w0, w) and 
x' = (w'w , «/). We now expand the Fourier basis into spherical harmonics, using 

e** = An V i'j,(qw)Y; m (0 x )Y lm q ) , (134) 



lm 



where X and q are the direction angles of * and q, respectively. The integration over the angle 
q can then be carried out by means of the orthogonality relation of the spherical harmonics, 



and ( |133[ ) turns into 

(ifr(0)ifrm) = - Zj J dw J dw 



an 



X qZdqP^MqwMqw^Y^YUtf). (135) 
Jo 

Expanding the effective lensing potential into spherical harmonics, we can write 

m = Yj ' ^ = f <?emY; m (0) d36) 
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and define the angular power spectrum by 

(tylm^Vm) = <5u>S mm 'Cf . (137) 



With ( |135[ ), this gives 

Cf = - d W du/^ ! — r q 2 dqP^q)j l (qw)j l (qw'). (138) 

T JO Jo K> s l«' Jo 

This expression allows a further simplification if the typical scales in the large-scale 
structure, 2n/q, are much smaller than the cosmological distances, w and w'. Then, qw :» 2n, 
and the spherical Bessel functions ji(qw) and ji(qw') vary very quickly compared to P^iq). 



We can then use 

r*ca 

q 2 dq ji(qw)ji(qw') = ^S(w - w) , (139) 
jo 2w 2 

evaluate P^iq) at q = l/w where the spherical Bessel functions ji(qw) peak, and approximate 



Jo 



1= f < 

Jo 



w s — w\ _ / / 



cr= | da,(-^| /> (-| . (140) 



This expresses Limber's H2761 127711 approximation. 

Since the Newtonian potential is given by the density contrast 6 in terms of the Poisson 
equation 

3 

V 2 = = ?"m05 , (141) 

2a 

the power spectra of the potential and of the density contrast are related by 

W = 4^1t) (142) 
and the power spectrum of the effective lensing potential becomes 

For the convergence k, we use did 1 Y[ m (6) = 1(1 + l)Yi m (ff) in Eq. ( 136) and find, for / s> 1, 

c;- 9 -^) 4 ^ J*Z!§ p,(L) . ( ,44) 

4 V c I Jo \w s a(w)j \w) 

As we shall see, this is identical to the power spectrum of the shear y, C] = Cf. Before we 
proceed, we need a digression to clarify the treatment of spin-weighed fields on the sphere. 



3.2. Basis vectors and spin-s fields on the sphere 

So far, when we had to introduce coordinates, we used the flat-sky approximation and kept the 
treatment restricted to local domains in which ordinary, Cartesian coordinates could be used. 
This is often adequate, but it is instructive and sometimes necessary to extend the domain such 
that the curvature of the sky can no longer be neglected. 
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There are two bases commonly used on the sphere. One is the orthonormal coordinate 

basis 

ex=d e , e 2 = -^- (145) 
sm9 

with the dual basis 

e l =dO, 2 = sin(9d0. (146) 
In this basis, the metric has the components g = diag(l, 1). Often, the helicity basis 

e± = 4= Oi + ie 2) (147) 
and its dual basis 

0* = (0 1 + i# 2 ) (148) 

are more convenient because the basis vectors are eigenvectors under rotations 7?(</?) on the 
sphere, 

<R(ip)e ± = e+ iv e ± . (149) 
Since the metric has the peculiar form 

9tJ = (° l = # d50) 

in this basis, pulling the indices + up or down changes their signs. The only non-vanishing 
connection forms are 

0) + ,=-ur = - C ^r{6 + -e~) , (151) 
V2 

from which covariant derivatives are easily constructed. 

Because its basis vectors are eigenvectors under rotations, the helicity basis is particularly 
useful for representations of spin- weighed quantities. A function s f is said to have spin s if 

J -» e" lv V (152) 

under right-handed rotations of the coordinate frame by an angle <p. Let T be a tensor field 
of rank s defined everywhere on S 2 . It defines a spin- s field on S 2 when applied to e_, and a 
spin-(-5) field when applied to e+, 

s t = T(e-,...,eJ), - s t = T(e + ,...,e + ) . (153) 

The covariant derivative V7~ of T is a rank-(^ +1) tensor field which defines a spin-(^ +1) 
field 

6 s t = -V2(V7~)(e-,...,e-,e_) (154) 

and a spin-(-5 - 1) field 

tf.,t = - V2(VT)(e + , ...,e + , e + ) . (155) 

The operator 5 is called "edth" B65H . It raises the spin by unity, while its adjoint & lowers 
the spin by unity. 
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Spin-s fields on the sphere can be decomposed into the spin- weighed spherical harmonics 
s Yi m defined by 



^<M=J ( ' ±m r +i) „ t Y,„, (156) 



such that 



±Ylm = 2 V(T^y. T2Y,m ' v+v - 7/m = Ylm = v - v+F/m ■ (157) 

The covariant derivatives V+ are short-hand notations for V e± . 
3.3. E and B modes 

Let now P be a symmetric, trace-free, rank-2 tensor field on S 2 . In the coordinate basis {e\, e 2 ), 
V may have the components 

This tensor field defines the spin-(+2) fields 

P(e ± ,e ± ) = {Pijff <g> 9 j ) (e ± ,e ± ) = \{Q + iU) =: T2 p . (159) 



The shear tensor Y in the coordinate basis {e\,e2\ introduced in Eq. (58 ) is an example for such 
a tensor field, with j\ = Q and y 2 = U. Another example is the linear polarisation tensor, 
where Q and U represent the usual Stokes parameters. 

The fields ±2 p are thus conveniently decomposed into spin-(+2) spherical harmonics, 

±2P = 1 -(Q + W) = Yj Vmm ±%Yi m , (160) 

Im 

with expansion coefficients p ±2 j m . Applying V 2 to - 2 p will raise the spin from -2 to zero, 



2 1 1(1 + 2)1 

V_- 2j p =: q- = - 2^ P-2M - J (i_ 2 y y ''« ( 161 ) 



Im 

72 



according to Eq. (157 ), while applying V+ to 2 p will lower the spin from 2 to zero, 



V +2 p =: q + = - 2_j P2 ' ,m ^ (j_2)\ lm ' ( ^ 

Since the two fields q± both have spin-0, they are independent of the orientation of the 
coordinate frame in which they are measured. Now we have, from Eqs. ( |161 ) and (162), 



f * 1 

]±,im = I <iQ.q+Yi m = -p± 2 j m 



and thus 



P±2,/m = 2 A / 2 ^ g ± ,; m . (164) 
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Imagine now observing the celestial sphere once inside-out and once outside-in. North 
will remain North, but East and West will be interchanged by this change of viewpoint, i.e. the 
coordinate basis {e\,e 2 } will transform into {e\, -e 2 }. Since 

Q = P(e 1 ,e 1 ), U = P(e u e 2 ), (165) 



with the tensor Eq. ( |158[ ), 

Q^Q, U^-U (166) 

under this parity transformation. 

Since the tensor components Q and U can be combined from the spin-+2 fields ±2 p as 

Q = ( 2 p + -2P) , U = -i( 2 p- - 2 p) , (167) 
the linear combinations of spherical-harmonic coefficients 

^E,lm = ~ \P2M + P-2,lm) , &B,lm = _ i (Pljrn ~ P-2Jm) (168) 

define a parity-conserving £-mode and a parity-changing 5-mode. The terminology reminds 
of the electric and magnetic fields because the electric field, as a gradient field, is invariant 
under parity changes, while the magnetic field, as a curl field, is not. A simple example for 
a 5-mode is a cyclone on Earth: a satellite looking at it from above will see it with opposite 
parity as a sailor looking at it from below. 

With Eq. (|164[), we obtain the spherical-harmonic coefficients 



a t.hii - _ *' (/ + 2)! ^ + ' lm 1-' lm ' ' 



(/ - 2)! 

ciBjm = - 2i J - - — (q +Jm - q^ hn ) . (169) 

from Eq. ( 168[ ). They define the E and B modes 

E0) = Yj a E,imYi m (0) , B0) = ^] a BM Y lm (6) . (170) 

Im Im 

While this decomposition is completely general for spin-(+2) fields on S 2 , we know that 

±2j P = V 2 (A (171) 
for gravitational lensing. Decomposing the lensing potential into spherical harmonics, 

^ = 2_j ^lm Y lm ' ( 172 ) 



Im 



we find 



, 1 1(1 + 2)1 

= - 2_j \7[Z^i^ ^ Ylm ' (173) 

Im } 



and thus, by comparison with Eq. (|160[), 



1 (Z + 2)! 

p±2M = iijr2y* lm - (174) 
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Therefore, as,im = 0, showing that lensing alone cannot create a 5-mode distortion pattern. 

Simplifications are possible in the flat-sky approximation, when we can set sin 6 = 1 in 
Eq. ( 145 ). Then, we introduce Cartesian coordinates #i = 6 and 62 = 4> an d have 

Vi2 = d 12 , V + = (d l + id 2 ) , 5 = fl 1 +ifl 2> tf = d-id2. (175) 

The decomposition into spherical harmonics can be replaced by Fourier transforms, 



such that 



f0) = YflmYlrn ~* f ^fti^ , 

V 2 J0) = - J ^L{l]-l\ + 2ihl 2 ) f{f)J 



(176) 



(177) 



for any field f(9) on the sphere. Identifying now / with Q + ill in Eq. ( 177 ) shows that 

q±M -> ~\ (/? ~lj)Q + hhU + i hhQ ll) U . (178) 

The flat-sky approximation is valid on small angular scales only, i.e. for /» 1. In this limit, 



a -2)1 ^ 1 

(l + 2)\ ~ I 2 
and Eqs. ( |164 ) and (169) give 



/?-Z 2 



Ah 



Ah 



1 2 - 1 2 

a EM -» -^fi - 2-^£/ , a BM ^^Q+ ±j±U 



Introducing the angle <p between I and 9, we thus find 

E - Q cos 2(p - U sin 2(p , B = Q sin 2(p + U cos 2^> 
for the Fourier-transformed E and 5 modes in the flat-sky approximation. 



(179) 



(180) 



(181) 



3.4. Shear correlation functions 

After this detour through spin- s fields on the sphere, we can now return to the treatment of 
shear correlation functions on the sky. Clearly, the shear components depend on the local 
coordinate frame and must be parallel-transported along well-defined curves when they are to 
be compared at two different locations on the sky. To avoid complications and ambiguities, 
it is much more appropriate to first combine them into spin-0 fields which are independent of 
the local coordinate orientations, and then compute their correlation properties. The preceding 
discussion tells us how to do this in a systematic way. 

Locally, we can return to the flat-sky approximation and combine the shear components 
712 into the spin-+2 quantities 

7 = 7i+m, 7* = 7i - i72 ■ (182) 
Expressed by the edth operator in the flat-sky approximation 5 = d\ + id 2 , 

K=~tf5ifr, y=-6V, F = , G = ^5>, (183) 
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which shows immediately that the F and G flexions as defined in Eqs. (69 ) and (70) are spin-1 
and spin-3 fields, respectively. 

The (negative) determinant of the shear tensor, 

Irl 2 , (184) 

has spin-0 and is thus independent of the orientation of the local coordinate frame. In the 
flat- sky approximation, we can Fourier transform the shear and find 

Z 4 



- det T = y j + y\ = 77* 



77 = 7i+7 2 = 4 W* 



KK 



(185) 



which shows immediately that the shear power spectrum defined in this way is identical to the 
convergence power spectrum, 

Cj = {yf) = q , (186) 



which was already given in Eq. ( 144 ). Other shear power spectra can be defined by converting 
the shear into spin-0 fields in alternative ways. The E and B modes defined in Eq. ( 170[ ) and 
expressed in the flat-sky approximation in Eq. ( 181 ) have spin-0 by construction. According 
to Eq. ( 1 8 their products are given by 

EE* = cos 2 2(pyi y \ + sin 2 2ipy 2 y* 2 and 
sin 2 2ipyiy\ + cos 2 2(py 2 y* 2 , 



BB* = 



(187) 



where 



7\7\ = 



(/? - 1\) 2 



(j/ijf* = cos 2 2(pkk* and 



Obviously, 



7272 = (hh) 2 ftft* = sin 2 2ipkk* . 



EE* + BB* = kk* = yf and EE* - BB* = cos 4(pyy* 



(188) 



(189) 



The two-point correlation functions of convergence and shear are given by the inverse Fourier 
transforms 

Idl 
— ( 
2n 

or the redundant definitions 

Idl 



-qum = &(0) 



(190) 



U0) = 



( c a 

Jo ~^ C * hm { + ) 



Idl , 

—quw) (-) 



(191) 



Examples for these shear correlation functions are shown in Fig. [5] 
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Figure 5. Shear correlation functions £ ± (9) as functions of the angular separation in arc 
minutes. The standard cosmological parameters given in Tab. [T] are assumed, and the source 
redshift is z s = 1. 



3.5. Gravitational lensing of the CMB temperature 

Gravitational lensing of the CMB has the effect that the relative temperature fluctuation 
observed in direction 9, 



T{9) - T 



(192) 



is the intrinsic temperature fluctuation in the deflected direction 9 - a, 

T oh M = T$-c$0)] . (193) 
Since typical scales in r are substantially larger than those in a, we can Taylor-expand 



Eq. ( 193 ) and write 



Tobs(^) = t0) - a ■ Vt(9) + l -4r^^j ■ 



(194) 



2 dOidOj 

These terms are more easily expressed in Fourier space. Since the deflection angle is the 



gradient of the lensing potential Eq. (51 ), 



(195) 



in terms of the Fourier transforms \j/ and t. Consequently, the scalar product a ■ Vr can be 
written as 

< 2 ' r ^ .... 

(196) 
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which has the Fourier transform 
d 2 /. 



r d% 

J (2nY 



Likewise, we can write 

12 



h-a-hmm-h) . (197) 



1 <9 2 T If d 2 /! 

i^hJm'** 1 (198) 



whose Fourier transform is 
12 



1 f tt4 f 44$ • • (**- £ - ^)]t(fi)^(4)<A(f- h - 12) . (199) 

2 J (2tt) 2 J (2^) 2 



Imagine now we Fourier transform Eq. (194) using the expressions Eq. (197) and 



Eq. (199), square the result, keep terms up to second order in the lensing potential, and 
average them over the sky to obtain the power spectrum of the observed CMB temperature 
fluctuations. In increasing order of powers in the deflection angle, the first term squared yields 
the intrinsic CMB power spectrum. The second term squared, and the product of the first and 
the third terms, lead to a convolution of the power spectra for r and iff. The product of the first 
and the second terms vanishes because the mean deflection angle does, and other terms are of 
higher than second order in iff. A somewhat lengthy, but straightforward calculation gives 

12; 



where 



q obs = c[ (i - 1%) + / iA. p . ( f_ /of c{c^ , (200) 
i r°° 

R f = — PdlCf * 3 x 10~ 7 (201) 
3^r Jo 

quantifies the total power in the deflection angle. 

We thus see two effects of lensing on the CMB temperature power spectrum 043811 . First, 
its amplitude is lowered by a factor increasing quadratically with the multipole order. Second, 
it is convolved with the power spectrum of the lensing potential, which gives rise to a slight 
smoothing. For / > 2000, the factor (1 - l 2 R^) becomes negative, which signals the break- 



down of the Taylor approximation (194). There, however, the amplitude of CT is already so 



low that this is practically irrelevant. For large /, we can thus ignore the first term in (200) and 
approximate (/ - l{) ~ I in the second term. This gives 



: (202) 



i r°° 

q obs ~i 2 cfR T with r t = —\ Pdiq « io~ Vk 

37T Jo 

Lensing thus creates power in the Silk damping tail of the CMB [ 3361 . 
3.6. Gravitational lensing of the CMB polarisation 

Since Thomson scattering does not produce circular polarisation, the polarisation of the CMB 
can be described by a symmetric, trace-free, rank-2 tensor field V whose components in the 



coordinate basis are the usual Stokes parameters Q and U as given in Eq. (28 ). As described 
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before, and entirely analogous to the shear, we thus decompose the polarisation into E and B 
modes which, in the flat-sky approximation, are given by Eq. ( [53] ), which can be combined 
into the form 

S + i!B = e ±2i *(Q±iU) (203) 

To avoid confusion with the E and B modes of the gravitational shear, we use curly symbols 
for the fi and S modes of the polarisation. Since no S-mode polarisation is expected to have 
been produced when the CMB was released, we can write 

Q + iU = e+^fi . (204) 

The following steps are conceptually simple, but require lengthy calculations. First, we 
repl ace t he temperature co ntrast r by the combinations Q + ill of Stokes parameters in 
Eq. (194). According to Eq. (204), this implies replacing f(Zi) by e + m S(Ji) in the expressions 



Eq. (197) and Eq. (199). The results need to be multiplied by e ±2l< ^ to find the observed, 
i.e. lensed, combinations (£ + iS) obs . Having obtained those, we can form the power spectra 



((£ + i£)(S + iff)*) 



Cf + Cf 



and 



(205) 



(206) 



((£ + rS)(£ - iS)*) = Cf - Cf . 
Their sum and difference yield the separate power spectra Cf and Cf , 

+ / [ fl(f " ^ C ° S2 2(01 " 0)C ^I C ' ' 

C is = / H [Ul- h)f sin 2 2(0! - <p)Cl f Cf { . 

Thus, the observable fi-mode polarisation power spectrum differs from the temperature power 
spectrum Eq. ( 200 ) merely by the phase factor cos 2 2{(p\ -0) in the integral, where (p and (f>\ are 



(207) 



the angles enclosed by I and l\ with the x\ axis. This phase factor differs from unity when the 
Fourier modes of the lensing potential and the intrinsic fi-mode polarisation are not aligned. It 
is very important that any such misalignment will create S-mode from fi-mode polarisation, 



as Eq. (207 ) demonstrates: when <p\ ^ <p, the phase factor sin 2 2{<p\ - <p) differs from zero, and 
the observed 5-mode power spectrum is different from zero although any intrinsic S-modes 
were explicitly ignored. Thus, we have arrived at a third effect of gravitational lensing on the 
CMB: lensing will create S- from fi-mode polarisation [151511 . 

Structures in the CMB are typically much larger than structures in the gravitationally- 
lensing density field. Thus, we can take the limit / l\ and approximate 

Cf obs « / sin2 2(0i - <PKCI = J C\CI , (208) 

which becomes independent of /. 
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3.7. Recovery of the lensing potential 

We do not know the appearance of the unlensed CMB, but lensing produces characteristic 
patterns on the CMB which may reveal its presence. In particular, lensing mixes Fourier 
modes and thus correlates the formerly unrelated modes 

t obs (Z) and f obs (7 + L). (209) 

To lowest order in the lensing potential iff, 



T Ob s(0 = t(/)- f ^m-tmim-h), (210 

J (2tt) 2 



so the expected correlation of different modes turns out to be 

(tobs(/)r; bs (/*- fy = (2n) 2 CJ6(t) (211) 



+ 



L(L-T)C r u +L-lC} 



To lowest order in (jj, lensing creates off-diagonal terms in the CMB temperature power 
spectrum, i.e. it couples different fluctuation modes which would otherwise be independent. 
This motivates the construction of a filter w(l, L), of which we require that the convolution 

N(L) J ^[f obs (/1f; bs a-L)]»;(Z:L) = (Aest(L) (212) 

return an unbiased estimate fy e $i(L) of the lensing potential. The requirement that the filter be 
unbiased, 

(<Aest) = <A , (213) 
puts a normalisation constraint on the prefactor N(L), 

L(L-l)C^ + L-iq 



, -> r d 2 i 



w(l,L). (214) 



The shape of the filter w(l,L) follows from the requirement that the significance of its 
detections be optimised, i.e. that the variance 

((«A est - <A) 2 ) = (|«A est | 2 ) + (|«AI 2 ) - (|<Aestl 2 ) (215) 
be minimised. Inserting the expression Eq. (212) yields 

(fcWestOL')) = N 2 (L)6(L-L') J ^La h C r ^_ n w%X) . (216) 

Varying this with respect to the filter w{l,L) under the bias constraint Eq. (213), taking into 
account that the normalisation factor N(L) also depends on w(f,£), and setting the variation to 
zero returns the filter 

L{L-l)0 + L.IC] 
w(l,L) = , (217) 

\ti-t\ ' 

whose application to the squared temperature field will return an unbiased, optimised estimate 
of the lensing potential H212H . This demonstrates one example for the recovery of the lensing 
effects on the CMB (see also AM 15751 ). 

These remarks conclude the formal part of the review. We shall now proceed discussing 
applications of gravitational lensing to cosmic mass distributions of increasing scale. 
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4. Structure and contents of galaxies 

4.1. MACHOs and planets 

We know that the most of the matter in the universe is dark, but we do not know what this dark 
matter consists of. We know that it must not interact electromagnetically, because otherwise 
the cosmic microwave background would show temperature fluctuations on the level of 1(T 3 K 
rather than 10 5 K. We also know that the dark matter must be cold in the sense that the 
velocity of its constituents must be small compared to the speed of light, because otherwise the 
large-scale distribution of the galaxies would be different. It is likely that this cold dark matter 
is composed of weakly-interacting elementary particles, but it could equally well consist of 
compact objects like, e.g. low-mass black holes. 

Gravitational lensing provides one way to test this possibility. Our galaxy, the Milky 
Way, is expected to be embedded into a halo which also predominantly consists of dark matter. 
If that dark matter was composed of compact objects rather than elementary particles, lines- 
of-sight through the Galaxy would occasionally pass close to one of those. They would act 
as point-mass lenses on sources in their background H112] . Although their image splitting 
would be substantially below the detection threshold, they would cause a well-measurable 
magnification 0761 . 

Quite independent of the mass spectrum of these hypothetical compact objects (called 
MACHOs^]), the probability of any one of them causing a microlensing event at any instant 
of time is of order (v/c) 2 , where v - 220 km s -1 Mpc -1 is a typical velocity for the stars 
in the Galaxy. Consequently, this microlensing optical depth is of order 10~ 6 . Finding its 
magnification signature thus requires of order 10 6 light curves to be monitored. Originally 
perceived more as science fiction, projects were carried out which observed sufficient numbers 
of stars in the Large and Small Magellanic Clouds (LMC and SMC, respectively) with 
sufficient accuracy and time sampling for detecting some microlensing events among the 
overwhelming signal from variable stars. 

Analysing data taken from 11.9 million stars over 5.7 years, the MACHO project [OQ 
found 13-17 events, while between 2 and 4 were expected from known stellar populations in 
the Milky Way and the LMC. The microlensing optical depth deduced from lensing events 
lasting between 2 and 400 days is t = 1.2^3 x 10" 7 . This implies that between 8% and 
50% of the Milky Way's halo can be composed of MACHOs (at 95% confidence), whose 
most likely mass ranges between 0.15 and 0.9 M Q (cf. Fig. [6]). Consistently, the EROS project 
H267B found, based on observations of LMC and SMC, that MACHOs cannot dominate the 
Galactic halo if their masses are < 1 M©. They find that the halo mass fraction in MACHOs 
is < 20% for MACHO masses between 10 7 M Q and 0.1 M Q (at 95% confidence). 

Thus, although MACHOs have been detected between us and the Magellanic Clouds, 
they are insufficient for explaining all of the Milky Way's dark mass. These MACHOs can 
in principle be anywhere between the source stars and the observer, i.e. in the dark halos of 
the Milky Way or of the Magellanic Clouds lfTMl4T3ll5M[T55l[TT6ll . It had been speculated 

§ acronym for massive compact halo objects 
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Figure 6. Example of likelihood contours obtained from the MACHO experiment for 
one specific model for the Milky Way halo, with different cuts applied to the candidate 
microlensing events (left and right columns). The abscissa is the fraction of the halo mass 
contained in MACHOs, the ordinate is the MACHO mass. The contours show the 60%, 90%, 
95%, and 99% confidence levels, (from Q) 

that self-lensing of stars in the LMC by its own stars might suffice for explaining the observed 
optical depth [7], but later studies showed that certainly not all of the LMC lensing events can 
be explained as being due to stars in the LMC. Rather, the LMC needs to be embedded into an 
extended halo HI 65112301 12971 . It seemed thus safe to conclude that microlensing experiments 
have confirmed that the Galaxy and the Magellanic Clouds must have extended dark haloes, 
only a fraction of which can be composed of compact objects of stellar and sub-stellar mass. 

With microlensing monitoring programs targeting the Andromeda galaxy M 31, the 
situation became less clear. Several early detections of microlensing candidates H3811 
14051 1 100L |480] proved the feasibility of such surveys and were tentatively interpreted as 
demonstrating microlensing in the halo of M 31. Detailed analyses of meanwhile more than 
a dozen candidates arrive at different conclusions regarding the contribution by self-lensing. 
While H219II argued that self-lensing should only occur near the centre of M 31, H101II found 
that all candidates are consistent with self-lensing and rejected the hypothesis at the 95 % 
level that > 30 % of the halo mass could be contributed by MACHOs. Conversely, [1661 had 
used six short-duration events to conclude that the MACHO contribution should be > 20 % if 
their mass was around 1 M Q . The main differences between these analyses are the models for 
the distributions of mass and stars in M 31 and the selection of microlensing events. Taking 
finite source sizes into account makes self-lensing less probable 1140611 . Low variability in 
individual light curves of 29 multiply imaged quasars supports the impression that MACHOs 
contribute at most little to microlensing in galactic halos H319H . 

Even MACHOs in distant lensing galaxies have been constrained by means of 
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microlensing. If a quasar is multiply imaged by a galaxy, stars and hypothetical MACHOs 
in that galaxy can act as microlenses on all images individually, giving rise to independent 
brightness changes in them. From the amplitude and the time scales of these brightness 
variations, assuming a plausible model for the size of the quasar source, upper limits can 
be placed on the fraction of the galaxy mass in compact objects within a given mass range 
H495L Turning this argument around, the size and the structure of the quasar engine can be 
constrained assuming microlensing by the galactic stars Il4l lllli 




days since 31.0 July 2005 UT 

Figure 7. Microlensing lightcurve showing the signature of a planet with 5.5 Earth masses, 
(from E3) 

A fascinating application of microlensing was first discussed by [130411 who showed that 
binary lens systems may leave detectable signatures in microlensing lightcurves even if one of 
the companions is a planet. This is strongly favoured if the planet's distance from its host star 
approximately coincides with the star's Einstein radius [156|. While the first planets detected 
in this way had Jupiter-like masses [|47ll479ll . even Earth-like planets have been found ( OBI , 
see Fig. [7]). To date, ten planets in nine planetary systems have been detected by means of 
microlensin^H} 



4.2. Galaxy density profiles 

Well over 100 cases of strong lensing by galaxies are now known. Most of them have two 
or four images, but a few have higher image numbers. Image splittings, typically of order an 
arc second, allow the projected lens mass to be constrained which is enclosed by the images. 
However, it turns out to be surprisingly difficult to constrain the mass profile from strong 
gravitational lensing alone. Essentially, multiple images constrain the average surface-mass 
density in an annulus bounded by the images. 

Using many lens systems, and assuming their mass profiles to be self-similar, it becomes 
possible to trace the average surface mass density at different radii, and thus to map out 



see the data base at http://exoplanet.eu/ 
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the density profile. Analysing 22 galaxy lenses, and adapting a mass model composed of 
a concentrated component representing the light and a power-law component representing the 
dark matter, [140911 found that the slope of the density profile is very nearly isothermal, with 
a double-logarithmic slope of n = 2.06 + 0.17 (isothermal has n = 2). Models in which the 
mass traces the light and is therefore more centrally concentrated fail at the 99% confidence 
level H410L If the dark matter follows the NFW density profile, (22 + 10)% of the matter 
inside two effective radii has to be dark. A weak trend is also seen in the mass-to-light ratio, 
M/L oc L 014 -o i2, consistent with the fundamental plane of elliptical galaxies. H478I agree that 
the density profiles are nearly isothermal, but find a somewhat larger scatter. They confirm 
that lensing galaxies in which light traces mass are ruled out, and find a dark-matter fraction of 
between 15% and 65% within the effective radius. Combining stellar kinematics with lensing 
in 15 lensing galaxies of the Sloan Lens ACS Survey (SLACS, 114611 ), [14770 find isothermal 
mass distributions within the Einstein radii. The density profiles found do not depend on 
environment H476L 

The picture becomes more complete when constraints from strong and weak lensing are 
combined. While multiple images produced by strongly lensing galaxies support isothermal 
mass distributions when interpreted with power-law density profiles, they probe the density 
profiles only in the narrow range of radii enclosing the multiple images. Weak lensing can be 
observed to much larger distances from galaxy centres. The combined analysis of 22 elliptical 
lensing galaxies by fll44] is one example further discussed below. Weak lensing was used to 
improve the mass models for the double quasar 0957+561 H130H and to refine estimates of the 
Hubble constant from its time delay H3571I1181 . EH studied the environment of eight lensed 
quasars by means of weak lensing and found indications for nearby galaxy groups in five of 
them. 

We have seen in the introduction that gravitational lenses are expected to produce an odd 
number of images. In contrast, all but very few observed galaxy-lens systems have an even 
image number, most of them either two or four. The missing images are expected to be faint 
if the central density profile of the lensing galaxies is steep enough, thus their absence can 
be used for constraining the central concentration of the lensing mass distributions. Based 
on this argument, [141 111 find that inner mass distributions of lensing galaxies cannot be much 
shallower than isothermal. Conversely, H504II use a lens system in which a faint, central image 
has been found to constrain the mass of the central black hole in the lensing galaxy to be 
< 2 x 10 s h~ l M e . From the general absence of faint, central images, 1124111 concludes that 
the central mass profiles of lensing galaxies must be more concentrated than CDM alone 
predicts. Central black holes may reconcile CDM density profiles with even image numbers 
only if they are about an order of magnitude more massive than expected from the relation 
between black-hole and bulge masses. 

Time-delay measurements in multiple-image systems promise constraints on the Hubble 
constant, provided a sufficiently accurate mass model for the lens is known. Conversely, 
considering the Hubble constant as known, time-delay measurements can be used as further 
constraints on the lensing density profile. Values for the Hubble constant derived this way 
tended to be lower than those, e.g. obtained from the HST Key Project (e.g. H123II ). but lens 
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models constructed upon a large number of constraints yield values which are very well 
in agreement with other determinations. For instance, (126311 find H = 75^kms _1 Mpc -1 
from time delays measured in the quadruply lensed quasar B 1608+656. From detailed 
modelling of the same lens system, H4671 derive Ho = 70.6 + 3.1 kms" 1 Mpc -1 for a 
standard cosmological model, and H = 69.0^^ km s" 1 Mpc -1 when combined with the 
WMAP constraints. In the latter case, the equation-of-state parameter is constrained to be 
w = -0.94+Q-Jg. Non-parametric mass models for lensing galaxies provide an interesting 
approach alternative to the parameterised density profiles. With this technique, H413I 
construct non-parametric mass distributions for 10 galaxy-lens systems and derive H = 
72.8^ 3 ° 5 kms~ 1 Mpc" 1 , while 0781 find# = 66^kms _1 Mpc" 1 from 18 lens systems. 

An interesting problem with interpreting time delays in four galaxy-lens systems was 
pointed out by [260]. It turns out that the time delays between multiple images is essentially 
determined by the mean surface-mass density in an annulus around the lens centre bounded 
by the images [259]. An estimate for that mass density can also be obtained converting the 
visible light to mass, assuming typical values for the fraction f b of matter that condenses 
into stars. Adopting ft, ^ 0.02 in accordance with local observations works well with near- 
isothermal mass models, but yields substantially too low values for the Hubble constant, 
Hq = (48 + 5)kms _1 Mpc" 1 . Conversely, values for the Hubble constant agreeing with the 
HST Key Project result, H = (72 + 8) km s~ ! Mpc -1 , are compatible with the measured time 
delays only if lens models with constant mass-to-light ratios are adopted, which are otherwise 
ruled out H260H . Using 14 time delays and a prior for the Hubble constant derived from other 
cosmological observations, f 10511 find that the density profiles of the lensing galaxies can only 
barely be compatible with being isothermal (see Fig. [8]). In a detailed analysis including 
various deviations from idealised lens models, [136811 finds H = (68 + 10) km s" 1 Mpc -1 
including systematic errors. 

There is thus the apparent problem that multiple-image configurations in galaxy lenses 
favour isothermal density profiles, while values for the Hubble constant derived from time 
delays are only in agreement with independent measurements if the density profiles are steeper 
than isothermal. A possible solution was pointed out by 113971 who showed that measurements 
of time delays and the Hubble constant can be reconciled if perturbations of the density 
profiles around isothermality are allowed. 

Altogether, the picture emerges that galaxy mass distributions approach isothermal or 
slightly steeper density profiles where baryonic physics dominates, and turn into the NFW 
profiles expected from numerical simulations for stable dark-matter dominated structures 
H3001 11441 130 IB (see Fig. [9]). Selecting galaxies by strong lensing may lead to substantially 
biased samples, though H302L while the spectroscopic selection defining SLACS Il46l 
identifies an unbiased galaxy sample [[81 I159II . 

4.3. Substructure in lensing galaxies 

Interestingly, the observed lenses with four images (quadruples) are about as abundant as 
such with two images (doubles), while they should only contribute 25% to 30% of the galaxy 
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Figure 8. Likelihood contours in the plane spanned by the Hubble constant and the double- 
logarithmic slope of the density profile assumed for the lens models, obtained from fitting 14 
time-delay measurements, (from |105|) 
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Figure 9. Surface-density profiles obtained from 22 strongly-lensing elliptical galaxies, 
compared to expected stellar and dark-matter components, (from [ 144|) 



lenses. The fraction of quadruples can be enhanced by satellite galaxies orbiting the main 
lens galaxies Il871l or by matter in their larger-scale environments H245L although the latter 
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explanation is potentially problematic because it also tends to lower inferred values of the 
Hubble constant. 

Axially- symmetric lens models are insufficient for modelling observed multiple-image 
systems. At least elliptical lens models are necessary. Embedding the lenses into additional 
external shear fields helps fitting observed image configurations, but typically more shear is 
required (10%- 15%) than the average large-scale structure can provide (l%-3%, H243I0 . This 
hints at the presence of asymmetries and substructures in the lensing galaxies. 
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Figure 10. MERLIN radio image of the multiply-imaged quasar B 1422+231 (courtesy of 
the JVAS/CLASS team, http://www.jb.man.ac.uk/research/gravlens/lensarch/). If image by a 
single, unperturbed lens, the two outer of the three bright images should together be as bright 
as the one between them. 



It is an interesting problem which caused much recent discussion that lens models are 
typically very successful in reproducing image positions, but fail in a large fraction of lens 
systems to explain the flux ratios between different images. A particularly obstinate and 
well-known case is B 1422+231 H264[ 120811 . see Fig. 10 This is most striking in situations 
where the source falls just inside a cusp of the caustic curve, in which case the sum of the 
signed magnifications of the three related images should vanish exactly. This expectation is 
frequently violated in real lens systems. 
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[130511 first suggested that substructure in the lensing galaxy could account for these 
anomalous flux ratios. While microlensing by the stars in the lens would less affect radio 
than optical fluxes because of their larger size, lensing by larger-scale substructures would 
equally change radio and optical flux ratios. Alternatively, CDM galaxy halos should contain 
sub-halos which may also account for anomalous flux ratios. Il52ll found similar modelling 
problems for images produced by a simulated lensing galaxy as for B 1422+231. IT78TI 
discussed that the sub-halo population of CDM halos produces perturbations of the magnitude 
required for explaining anomalous flux ratios. [133711 estimated that > 5% of the lensing 
halo mass must be contained in substructures and argued that elliptically deformed power- 
law models embedded into external shear are insufficient for most lenses. Similarly, 11991 
concluded that substructure comprising 0.6% to 7% of the lens mass, with a median at 2%, 
was necessary for reproducing the observed anomalous flux ratios, in excellent agreement 
with CDM halo simulations. They also estimated that the sub-halos should have masses 
in the range (10 6 . . . 10 9 )/i _1 M Q . H334M found that halo substructures with masses within 
(10 5 . . . 10 7 ) h' 1 M Q may explain the curved radio jet in B 1 152+199. 

In contrast to these arguments, HI 171 explicitly constructed smoothly deformed lens 
models which could well reproduce image configurations and flux ratios for most lens systems 
and argued that substructure in the lensing galaxies and smoothly deformed lenses are both 
viable explanations for the anomalous flux ratios. Along the same line, [289] and H395II 
showed that disks in lens galaxies can alter image magnification ratios considerably, while 
TTl\ pointed out that halos projected onto the main lens galaxy may also cause the observed 
magnification perturbations. 

From a somewhat different perspective, H418I discussed that decomposing lensing 
galaxies into microlenses has the most prominent effect when part of the lensing mass remains 
smooth. Specifically, they showed that so-called saddle-point images can be substantially 
demagnified in presence of microlenses. H4191I added that the magnification distribution of 
the macro-images depends on the mass spectrum of the microlenses, in contrast to earlier 
expectations. 

H2611I investigated various alternative explanations for the anomalous flux ratios, 
such as absorption, scattering, scintillation, uncertainties in the macro-model, and stellar 
microlensing, and arrived again at the conclusion that halo substructures remain as the most 
likely reason. The parities of the three images formed near a cusp point in the caustic is an 
important argument to exclude other alternatives to lensing. llFTI verified that numerically 
simulated galaxies can produce anomalous flux ratios as observed and emphasised the 
importance of the demagnification of saddle -point images. Although it appears doubtless 
that CDM halos contain sufficient substructure for sufficiently perturbing image flux ratios, 
such sub-halos must also appear projected onto at least one of the images. The probability 
for that is low. [130311 find in numerical simulations that the probability of finding suitably 
massive sub-halos in front of macro-images is only < 0.5%. Newer simulations all agree 
that the expected level of substructures in lensing-galaxy haloes is insufficient for explaining 
the anomalous flux ratios, even if high mass resolution and baryonic physics are taken into 
account H2921 131 1291 115 131 . Halos unrelated to the lensing galaxies, but projected onto them, 
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can have a substantial effect on flux-ratio anomalies H3351 13451 . There is agreement that 
perturbing masses must be compact 11751 1452H . even though smooth perturbations appear in 
sufficiently flexible lens models [|414] . In some cases, detailed modelling based on all lens 
components visible in detailed observations reproduces the measured flux ratios well without 
any further substructure H3181I2941I . 

The situation thus remains interestingly confused. It seems likely that the combination of 
dark-matter substructures, the tidal field of the environment, mergers and random projections 
provides a way out. H244 | l suggest that time-delay perturbations may add further constraints 
on halo substructures. 

4.4. Lens statistics 

The abundance of galaxy lenses has often been used for constraining the cosmological 
constant A. While early studies typically found upper limits of A < 0.7 (e.g. H2581 112110 . 
more recent investigations find values which are better compatible with other determinations, 
(e.g. [1791 1731 ). finding spatially-flat model universes with low matter density (Q - 0-3) 
preferred. The reason for this change is that gradually more realistic galaxy luminosity 
functions were used for estimating the expected number of lenses, rather than error-prone 
extrapolations of local galaxy number densities towards high redshift 1124211 . 

The Sloan Digital Sky Survey has allowed the definition of a homogeneously selected 
quasar sample (121711 from which cosmological parameters were derived [369]. Assuming 
a spatially flat universe, a value of f2 A0 = 0.74^-Jg was derived for the cosmological 
constant, where statistical and systematic errors were combined in quadrature. Allowing 
a free equation-of- state parameter gave w = and a matter-density parameter of 

^mo = 0-26 + 0.08 when combined with independent cosmological constraints. 

Halos are expected to have a continuous mass spectrum in universes dominated by cold 
dark matter, which is described by mass functions such as those derived by [392], [145 IB and 
H229H . Thus, one would expect a continuous distribution of splitting angles between fractions 
of an arc second to several ten arc seconds. 1135911 investigated whether the observed image- 
splitting distribution was consistent with expectations from CDM. They found observation 
and theory agreed if selection effects were taken into account. H257II found that the splitting- 
angle distribution in CDM is grossly incompatible with microwave-background constraints in 
a model universe with high matter density and vanishing cosmological constant, but that both 
could be comfortably reconciled in a spatially-flat, low-density CDM model. 

Occasionally, therefore, lens systems should be detected with splitting angles of ten or 
more arc seconds. [389] interpreted the absence of wide-separation lenses in the CLASS 
survey as a being due to low central mass concentrations in group- and cluster-sized halos. 
It was perceived as a further confirmation of the CDM paradigm when a quadruply imaged 
quasar was detected in the Sloan Digital Sky Survey with a splitting angle of 14.62 arc seconds 
12181 . for which 1150211 derived a lens mass of (5± 1)X 10 13 /r 1 M Q within a radius of 100 h~ l kpc 
based on a non-parametric lens model. 1137011 noted that the triaxiality of CDM halos must be 
taken into account in probability and mass estimates for the formation of wide- separation lens 
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systems, a theme which is familiar from studies of strong lensing in galaxy clusters. 

Clearly, cosmological parameters from the statistics of strong gravitational lensing by 
galaxies are generally no longer competitive compared to those based on observations of the 
cosmic microwave background because the uncertainties in lens models and sample selection 
are considerable. It should be kept in mind, however, that the cosmic microwave background 
does not independently measure the Hubble constant, but the expansion rate during the time 
of recombination. Hence, independent measurements in particular of the Hubble constant are 
and remain most important. 

4.5. Weak lensing of galaxies by galaxies 

Less distant galaxies can act as weak gravitational lenses on more distant galaxies. Their shear 
imprints a feeble tangential distortion pattern on the images of background galaxies which 
appear projected close to them. This weak signal is superposed on the intrinsic ellipticities 
and irregularities of the background-galaxy images and thus requires statistical techniques 
for its extraction. Il55ll first discussed the principal features of this effect and searched for 
it in a sample of galaxies, in which they separated background from foreground galaxies 
according to their apparent brightness. They could already infer that the shear profile of 
brighter galaxies was compatible with an isothermal mass profile with a circular velocity of 
v c = (220 + 80) km s _1 . They also placed a lower limit r* > 100 h~ l kpc on the halo size of the 
lensing galaxies. (1 10211 searched for galaxy-galaxy lensing in the Hubble Deep Field (North) 
and found a mean velocity dispersion for the lensing halos of cr v = 185^° km s _1 , and a weak 
lower limit on the halo radius. 

H426B devised a maximum-likelihood technique for efficient analysis of galaxy-galaxy 
lensing data which specifically took the redshift distributions of foreground and background 
galaxies into account. They applied this technique to numerically simulated data and 
calibrated its performance. [136011 and [1 14611 developed methods for detecting the weak- 
lensing signal of galaxies embedded in galaxy clusters. Applying their technique to the cluster 
CI 0939+4713, [1 14711 detected the shear signal of individual massive cluster galaxies. More 
recently, H3611I compared the weak-lensing signal of early-type, L* galaxies in clusters and in 
the field and found evidence for the cluster galaxies to be truncated, with a truncation radius 
shrinking with the density of the environment. [206] combined weak-lensing data on galaxy 
halos to show that they are flattened. 

Recent wide-field surveys also triggered an exciting development of galaxy-galaxy 
lensing. H131H used the Commissioning Data of the SDSS to infer that the tangential shear 
profile is compatible with a power law with exponent between 0.7 and 1.1, i.e. close to 
isothermal. They found a best-fitting circular velocity of v c = (150... 190) km s" 1 and a 
lower limit to the physical halo radius of 260 h~ x kpc. From the Las Campanas Redshift 
Survey, 045 81 deduced an isothermal tangential shear profile within 200 h~ x kpc and a circular 
velocity of v c = (164 + 20)kms~ 1 for L t field galaxies. They found a virial mass for the 
dark halo of a typical L* galaxy of (2.7 + 0.6) x 10 11 h~ l M Q . [150311 used data taken with 
the UH8K camera at the Canada-France-Hawaii telescope to measure galaxy-galaxy lensing. 
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They also found tangential shear profiles compatible with an isothermal slope and a rotation 
velocity of v c = 238*?Z km s" 1 for L t galaxies. They concluded that the mass-to-light ratio of 
L t galaxies in the B band is M/L (121 + 28) h M Q /L Q . According to the analysis of 22 square 
degrees of CFHT Legacy Survey data by [137911 , L* galaxies at redshift 0.3 reside in haloes of 
(1.1+ 0.2) x 10 12 h' 1 M Q and have an i?-band mass-to-light ratio of 173 + 34 h M Q /L Q . 

Combining the galaxy-galaxy weak- lensing signal obtained from SDSS data with the 
Tully-Fisher and fundamental-plane relations for late- and early-type galaxies, respectively, 
H439H found that the galaxy velocity profile must drop substantially towards the virial radius, 
which indicates a steep dark-matter profile. H1631I compared theoretically motivated CDM 
halo models with SDSS data and constrained the halo properties of galaxies with luminosities 
> L t . They constrained the virial mass of an early-type L* halo to M 2 oo = (5 ... 10) x 
10 11 h~ x M Q , and somewhat less for late-type galaxies, depending on the colour. They found 
a gentle increase of the mass-to-light ratio with luminosity, with M/L ^ 11 h M Q /L Q for late- 
type and M/L ^ 45 h M Q /L Q for early-type L t galaxies. 

[448J studied the cross -correlation between galaxies and mass from the galaxy-galaxy 
lensing signal detected in SDSS data. The wide area covered by the survey allowed 
constraining the correlation function out 10 h~ l Mpc. They find a power law with a correlation 
length of r (5.4 + 0.7) h~ x Mpc and an exponent of 1.79 + 0.05. The bias parameter turns out 
to be approximately scale-independent (see also H204H ). while H203I find the bias parameters 
to be gently increasing from Mpc to larger scales. Comparisons with theoretical expectations 
for the galaxy distribution relative to the dark matter find overall good agreement H49611 , except 
that the simulated mass-to-light ratio is somewhat too high [15 1411 . Satellite galaxies orbiting 
the lensing galaxies could be physically aligned with their hosts and thus mimic a weak 
galaxy-galaxy lensing signal. H192H estimated this possible contamination and constrained 
it to less than 15% at the relevant scales. 

The availability of huge surveys with sufficient depth and image quality for weak- 
lensing studies has opened new applications also for galaxy-galaxy lensing. The statistical 
results obtained (e.g. by H3001 11441 13011 ) on the density profiles of galactic haloes were 
already described above. Exciting examples of more detailed studies enabled this way are the 
measurement of mean masses of halos hosting active galactic nuclei, which found that radio- 
loud AGN reside in host haloes which are typically « 20 times more massive than for radio- 
quiet AGN [299J, an d the combination of galaxy-galaxy lensing with galaxy correlations to 
specify the mean mass-to-light ratio of the galaxies [|65l . 

5. Galaxy clusters 

5.1. Strong lensing 

Strong lensing in galaxy clusters was first detected by 1146 IB and [284J. They found extended, 
arc-like images in the galaxy clusters A 370 and CI 2244. Several explanations were proposed 
for these objects, among them gravitational lensing of background galaxies H377H . which was 
confirmed when the redshift of the arc in A 370 was measured and found to be substantially 
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Figure 11. This HST image of the cluster SDSS J 1004+4112 14461 shows arcs, multiply 
imaged galaxies and a quadruply lensed quasar with a maximum image separation of 14.62 
arc seconds 12181 . 

It was quickly recognised that gravitational arcs provided important information on 
the structure of galaxy clusters. It was unclear at the time how the dark matter was 
distributed and whether the X-ray surface-brightness profiles, which typically show a flat 
core of - 200 h~ l kpc radius, were representative for the dark-matter profiles. Arcs were soon 
found to reveal the following about clusters: (1) Cluster mass distributions cannot typically 
be axially symmetric, because large counter-arcs would otherwise be expected HI 601 126511 . 
(2) The substantial amounts of dark matter in galaxy clusters cannot be attached to the 
galaxies because arcs would then have much smaller curvature radii H1711 1381 . Particularly 
striking were the detections of "straight arcs" in two clusters H3831 13151 139011 because they 
visually demonstrated the need for substantial concentrations of dark matter with very high 
mass-to-light ratio 1123911 . (3) Clusters need to have steep density profiles, because arcs 
would be substantially thicker otherwise (1 17011 . For clusters to be strong lenses, their central 
convergence k has to be close to unity, but for arcs to be thin, the convergence at their locations 
has to be around 0.5. From cluster centres to the arc radii of typically 10" . . . 30", the k profile 
must thus fall by approximately a factor of two. Cluster core radii, if they exist, must thus be 
substantially smaller than the X-ray core radii, which was also confirmed by the detection of 
"radial arcs" lfT33l 13401 132TI 

Arcs allow cluster masses to be easily estimated. It was soon discovered that the masses 
obtained this way are very close to mass estimates derived from the X-ray temperature and 
surface-brightness profile. This is not obvious because gravitational lensing is sensitive to 
the mass irrespective of its physical state, while the interpretation of X-ray data requires 
assumptions on symmetry and hydrostatic equilibrium of the gas with the gravitational 
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potential well, if not on isothermality of the intracluster gas. This overall agreement being 
reassuring, a systematic discrepancy was soon revealed in the sense that masses derived from 
strong lensing were typically higher by factors of ^ 2 ... 3 than X-ray masses H508L 13431 151 IB . 
There are many more recent examples. Many find substantially discrepant mass estimates 
based on X-ray and strong-lensing observations [13911 11491 13461 1 1 661 11101 , while good 
agreement is found in other clusters F412ll53l . 

[|29l used numerical simulations to show that X-ray mass estimates can be systematically 
lower in merging clusters because their X-ray gas is still cooler than expected from their total 
mass, which is already seen by the lensing effect. This seems to explain the mass discrepancy 
at least in some clusters (e.g. H456[ |3~740. Asymmetries and cluster substructures also play an 
important role. Due to their relatively larger shear, asymmetric and substructured clusters are 
more efficient lenses at a given mass. Mass estimates based on axially symmetric models are 
thus systematically too high irT3~lll73H . 

[0 distinguished clusters with and without cooling flows and found an appreciable mass 
discrepancy in clusters without, but good agreement of X-ray and lensing mass estimates in 
clusters with cooling-flow. This supports the concept that well-relaxed clusters which had 
sufficient unperturbed time to develop a cooling flow are well-described by simple, axially- 
symmetric models for lensing and the X-ray emission, while dynamically more active clusters 
tend to give discrepant mass estimates; this was confirmed by [15 1011 . [129611 noted that the mass 
discrepancy is reduced if cluster density profiles are steeper than inferred from the X-ray 
emission. It thus appears that mass discrepancies can commonly be traced back to dynamical 
activity in unrelaxed clusters (see also 059)), but at least part of the disagreement also 
occurs because of model restrictions which, if removed, generally lead to better agreement 

nsi[TEiE3a. 

5.2. Cluster mass profiles 

Assuming mass profiles with cores, tangential arcs require small core radii as described above, 
but radial arcs require the cores to be finite [268, 128211 . Numerical simulations of CDM halos, 
however, show that density profiles flatten towards the core, but do not develop flat cores 
H363[ 13641 . |[T5l showed that radial arcs can also be formed by halos with such "cuspy" 
density profiles, provided the central cusp is not too steep. 

In principle, the relative abundances and positions of radial compared to tangential arcs in 
clusters provide important constraints on the central density profile in clusters H34UI3511I372II . 
Radial arcs are still too rare for successfully exploiting this method. Being much closer 
to the cluster cores than tangential arcs, they are also more likely to be confused with or 
hidden by the light of the cluster galaxies. Following H341L [141611 compiled a sample of 
clusters containing radial and tangential arcs and added constraints on the central mass profile 
from velocity-dispersion measurements in the central cluster galaxies. They demonstrated 
that, assuming axially-symmetric mass distributions, central density profiles have to be 
substantially flatter than those found in CDM simulations. However, even small deviations 
from axial symmetry can invalidate this conclusion and establish agreement between these 
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observations and CDM density profiles lfT9l 132611 . 

Attempts at modelling arcs with isothermal mass distributions are typically remarkably 
successful (see [256J for an impressive early example). This is all the more surprising as 
numerical simulations consistently find density profiles which are flatter than isothermal 
within the scale radius and steeper outside. In a very detailed analysis, H141II find that an 
isothermal core profile for the cluster MS 2137 is preferred compared to the flatter NFW 
profile. H458H constrain the core density profile in A 383 using X-ray, weak-, and strong- 
lensing data and find it more peaked than the NFW profile, but argue that this may be due to the 
density profile of the cD galaxy. Similarly, H255B find in a combined weak- and strong-lensing 
analysis of CI 0024+1654 that an isothermal mass profile can be rejected, while the NFW 
profile fits the data well. (145311 show that the strong-lensing effects in two clusters A 370 and 
MS 2137 can be explained similarly well by isothermal and NFW density profiles, leading to 
substantial uncertainties in derived cluster properties and magnifications. Strong lensing alone 
constrains cluster density profiles only close to cluster centres, leaving considerable freedom 
in the mass models. Certainly, the innermost cluster density profiles can be significantly 
influenced and steepened by baryonic physics. 

It is a puzzling and potentially highly important problem that strong-lensing analyses of 
galaxy clusters in many cases find that NFW density profiles well reproduce the observed 
image configurations, but with concentration parameters that are substantially larger than 
expected from numerical simulations ll58ll88l[T78ll4811l59ll . This so-called overconcentration 
problem is currently much debated, in particular because other studies find concentration 
parameters in the expected range, sometimes in the same clusters HI 671 12791 . Due to selection 
biases and projection effects, strongly lensing clusters should be among the most concentrated 
clusters, but some clusters such as A 1689 seem to be extraordinarily concentrated. Whether 
this is a significant contradiction to the ACDM model remains to be clarified. 

5.3. Arc abundance and statistics 

The mean density profile of galaxy clusters can also be constrained statistically because 
the probability for a cluster to become a strong lens depends sensitively on the mass 
concentration in its core H512H . [134011 suggested that the core densities of strongly lensing 
clusters could be enhanced by projection of elongated clusters along the line-of- sight. PTfl 
used a numerically simulated galaxy cluster to show that asymmetric, substructured cluster 
models are significantly more efficient strong lenses than axially- symmetric mass distributions 
because of their enhanced tidal field. Averaging over a sample of simulated clusters, Il30l 
quantified that the cross sections for arc formation could be up to two orders of magnitude 
larger for asymmetric than for axially symmetric cluster models of the same mass. 

HI 741 confirmed that structured lenses help understanding the observational results of 
[268 1, who detected six arcs in a sample of 16 clusters selected for their high X-ray luminosity 
as measured by the EMSS satellite, but argued that even more concentrated mass profiles 
than those used by [30] are necessary for explaining them quantitatively. IfPTll used samples 
of numerically simulated clusters to estimate the total arc-formation probability in different 
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cosmological models. Comparing their results with the data from [268 J, they concluded that 
only their cluster sample taken from a simulation with low matter density (£2 = 0-3) and no 
cosmological constant could well reproduce the measured high arc abundance, but the other 
three models failed badly. In particular, a flat cosmological model with Qq = 0.3 and Q A = 0.7 
produced an order of magnitude less arcs than observed. 

This so-called arc-statistics problem was disputed based on calculations using analytic 
models for cluster lenses [|9T] 12401 . which failed to reproduce the strong dependence on the 
cosmological constant claimed by [17J. The possible influence of cluster galaxies on the 
arc-formation efficiency of cluster lenses was investigated by HI 321 and H329B , but found 
to be negligible. [135211 confirmed that axially-symmetric mass models adapted to the X- 
ray emission do not produce a sufficient number of arcs. They found that using NFW 
profiles for the dark-matter profile helped, but the profiles required too high masses, and 
proposed that substructured mass distributions could be the solution. [328J adapted elliptically 
distorted lenses with NFW mass profile (see also H152I0 to numerically simulated clusters and 
found the analytic models inadequate for quantitative arc statistics despite the asymmetry, 
demonstrating the importance of substructures. 

[137 IB studied the strong-lensing properties of triaxial (rather than ellipsoidal) halos and 
found that they may well explain the high arc abundance [268. 12821 . provided their central 
density slopes are steep enough, with a double-logarithmic slope near -1.5. [149311 simulated 
the magnification probability for light rays propagating through a section of the Universe 
and found that the abundance of high-magnification events depends strongly on the source 
redshift. They attributed this to the exponential mass function of massive halos, which leads 
to a steep increase with source redshift in the number of halos suitable for strong lensing. 
Identifying the probability for highly magnified light bundles on random patches of the sky 
with the probability for finding arcs in massive galaxy clusters, they suggested this result 
as the resolution for the arc statistics problem. [1981 used numerical cluster simulations to 
estimate arc cross sections and found reasonable agreement with the earlier results of IfPTll . 
but arrived at a higher expected arc abundance because they inserted a higher normalisation 
for the number density of both X-ray clusters and background sources. 

H501II noted that the arc radii in clusters depend only weakly on clusters mass and 
suggested that massive cD galaxies may be the reason. However, H327II studied the effect of 
cD galaxies on the overall arc abundance and found it insufficient to remove the arc statistics 
problem. If the cosmological constant is replaced by some form of dynamical dark energy, 
structures tend to form earlier during cosmic history. Since cluster core densities reflect the 
mean cosmic density at their formation time, clusters thus tend to be more concentrated in 
dark-energy compared to cosmological-constant models. [1201 estimated the effect of higher 
cluster concentrations on arc statistics by analytic means. They found that dark energy may in 
fact increase arc abundances noticeably, but again not sufficiently for solving the arc statistics 
problem. 

Galaxy clusters at high redshifts are found to be remarkably efficient lenses HI 501 15 161 
even though they should be by far not massive enough for producing large arcs. A particularly 
impressive example is the cluster RX J 105343 +5735 at z = 1.263 which contains a large arc 
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from a source at z = 2.577 H474H . In this respect, it is interesting that the strong-lensing 
efficiency of clusters can be increased substantially and on a short timescale during a major 
merger H475H . As a subcluster approaches a cluster, the tidal field is increased, leading to a 
first maximum of the cross section approximately when the two virial regions touch. The 
cross section then slightly decreases and approaches a second maximum when the separation 
of cluster and subcluster is minimal. A third peak corresponding to the first is formed when 
the subcluster leaves the virial region again after the merger. During that process, the arc 
cross section can change by an order of magnitude or more on a time scale of ^ 0.1 Gyr. It 
thus appears that strong lensing can be a transient phenomenon at least in some clusters which 
would otherwise be not massive or concentrated enough. The dependence of the main merger 
epoch on cosmic history would then establish an interesting link between high-redshift, strong 
cluster lenses and the cosmological framework model. 

So far, all effects studied, including baryon cooling in cluster cores H3931 1494H and 
line-of-sight projection effects H394L returned moderate enhancements of the expected 
arc abundance. The main problem is now that recent measurements converge on a low 
normalisation parameter <x 8 ~ 0.8 for the dark-matter power spectrum, which drastically 
lowers the expected arc abundance H2751 11261 . Thus, while the strong-lensing properties of 
individual clusters seem to be sufficiently understood H210II . the arc-statistics problem persists. 
It is certainly related to the overconcentration problem and the problem of large Einstein radii 
[|60l and indicates that there is something fundamental we do not understand in the population 
of galaxy clusters. 

Besides improved methods for predicting arc abundances H 12711 . selection biases in 
existing samples of strongly lensing clusters must be understood H 12511 and the sizes of 
observed arc samples need to be enlarged H179II before further progress will be made. 
Automatic search algorithms for arcs H27 1112 101 14321 l63l calibrated on simulated images [133 111 
and applied to upcoming wide-field surveys will improve the situation in the near future. 

5.4. Other applications of strong cluster lensing 

If a cluster produces arcs from multiple sources at different redshifts, the lensing mass 
distribution remains the same, but the geometrical lensing efficiency is different for the 
arcs. Since this depends on cosmological parameters, these can thus be purely geometrically 
constrained from multiple-arc systems [280i 11391 . [146211 applied this technique to multiple 
arcs in the cluster A 2218 and found that a universe with critical matter density and no 
cosmological constant is excluded at > 4 a confidence from this single cluster. In principle, 
the abundance of arcs in clusters and the geometry of multiple arc systems in individual 
clusters is sensitive to the possible time evolution of the dark energy, but in view of 
the uncertainties in cluster mass models, substantive conclusions will be hard to arrive at 
0301 13251 12901 . 

For alleviating potential problems e.g. with the abundance of satellite galaxies, it was 
proposed that the dark-matter particles might interact with each other in another way than 
through gravity. Such a self-interaction would act as a source of isotropic pressure and thus 
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symmetrise and smooth mass distributions H342L Strong gravitational lensing, being very 
sensitive to cluster asymmetries, places a tight limit on the interaction cross section. Using 
numerical simulations, H332H showed that the strong-lensing efficiency of galaxy clusters 
would abruptly disappear if the specific self-interaction cross section was > 0.1 cm 2 g _1 . 

Finally, lensing clusters are frequently being used as cosmic telescopes, magnifying 
distant sources above the limits for photometry or spectroscopy. To give a few examples, ll69l 
used the magnification by A 370 to detect CO lines in arc sources. [135] and (138211 identified 
sources at z = 4.04 lensed by A 2390, H2541I found an object with z — 7 lensed by A 2218, and 
(138411 claim to have detected an object with z = 10.0 magnified by A 1835. The magnification 
by galaxy clusters as cosmic telescopes has recently been used in systematic searches for 
galaxies at very high redshift H4641 14571 14041 l54l and for increasing the resolution in detailed 
studies of distant galaxies K468L 1469 1 . 

5.5. Weak cluster lensing 

Apart from the occasional spectacular strong-lensing effects, clusters imprint a coherent weak 
distortion pattern onto the many faint and distant galaxies in their background. Those galaxies 
reach number densities of ^ 40 per square arc minute in typical images taken with large 
ground-based telescopes. The virial region of a typical galaxy cluster thus covers of order 
10 3 galaxies. Due to their intrinsically irregular shapes, lensing-induced distortions cannot 
be inferred from individual galaxies. Averaging over a few galaxies, however, the intrinsic 
ellipticities should average to zero, leaving the ellipticity caused by the gravitational shear as 
the average signal. 

As shown in the introduction, shear and convergence are both related through the scalar 
lensing potential. Knowing the shear thus allows the scaled surface-mass density to be 
reconstructed. (123711 were the first to show that cluster convergence maps could be obtained 
by convolving the measured shear signal with a simple kernel, opening the way to systematic, 
parameter-free, two-dimensional cluster studies. Their technique was immediately applied 
to the cluster MS 1224, for which H1191I found a surprisingly high mass-to-light ratio of 
800 h Mq/Lq in solar units, about a factor of four times the typical cluster value. 

Weaknesses in the convolution algorithm by H237H . such as its limitation to formally 
infinite data fields and weak shear, were discussed and removed by, e.g. H43 3 1143411 . Another 
technique for recovering cluster mass maps based on a maximum-likelihood approach 
was proposed by iBTTl . later augmented with maximum-entropy regularisation K435K and 
further developed by H309K . Based on this approach, algorithms have been developed for 
combining constraints from weak and strong gravitational lensing into unique, non-parametric 
reconstructions of cluster mass distributions Ii50ll4^1l64ll333ll . 

An algorithm for measuring the weak shear signal from data fields was first described 
and implemented by K238K . The Shear Testing Programme (STEP, H1861 13111 was launched 
to test and improve the accuracy of shear measurements from weakly distorted images of 
distant galaxies in the presence of several perturbing effects. 

These inversion techniques for cluster lenses have by now been applied to numerous 
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clusters, far too many to mention them all here. For most of them, the mass-to-light ratios 
turned out to be quite normal, i.e. M/L 250 . . . 300 in blue and M/L =t 150 .. . 200 in red 
colour bands, respectively. Some examples are It84l I200H and 0321. Statistical analyses of 
large samples of galaxy groups and clusters have allowed extending the mass range in which 
mass-to-light ratios can be probed. Values of M/L = 185 + 28 [380| in blue and M/L 250 
H278B in red spectral ranges have been obtained for galaxy groups with masses around 
(10 13 . . . 10 14 )M o . Cross-correlating the weak-lensing signal around a sample of 130,000 
brightest cluster galaxies from the Sloan digital Sky Survey with red light, H449L 145011 found an 
average mass-to-light ratio of M/L = (362 + 54) b 2 , where the bias factor b should be of order 
unity. Mass and light generally appear well correlated in weak-lensing clusters. The finding 
that mass followed light well in two of the three clusters in the A 901/902 super-cluster field 
while the third cluster showed a significant offset between the mass and the light [115 811 could 
not be confirmed by a follow-up study with HST data HI 841 . A very peculiar case is the cluster 
A 520, where the galaxies are seen displaced from the coincident dark-matter distribution and 
X-ray emission H295I . see Fig. [12 



Interesting phenomena appear in comparisons between the X-ray surface-brightness and 
the weak-lensing mass contours. While the X-ray emission follows the mass in many clusters 
(see H 14811831 [1991 [86l for examples), instructive deviations have been discovered. H293B find 
good agreement between surface-density and X-ray contours in the outer parts of A 2218, 
but deviations near the cluster centre, which they interpret as a sign of dynamical activity 
in the cluster. Several recent studies find the X-ray gas lagging behind the dark matter in 
merging clusters ll8Tl 13081 12281 . as expected for hot gas embedded into collision-less dark- 
matter halos. A particularly interesting case is the cluster IE 0657-558, called the bullet 
cluster, whose X-ray emission appears in between two galaxy concentrations and dark-matter 
distributions recovered from weak lensing Il82ll80l . see Fig. 13 Its appearance suggests the 
interpretation that two clusters have lost their gas by friction while passing each other in the 
course of a merger. Other such clusters have since been found [48J. 

If the dark-matter particles interacted with each other, such a separation between gas and 
dark matter would be suppressed. Thus, from gas lagging behind the dark matter in merging 
clusters, and from small dark-matter core radii, limits could be obtained for the self-interaction 
cross section of the dark-matter particles, typically finding values < (0.1 . . . 1) cm 2 g _1 
iH 13621 13081 . comparable to what 1133211 concluded from strong cluster lensing. Clusters 
like the bullet cluster have also been used to constrain the specific cross section for the self- 
interaction of dark-matter particles to cr/m < 1 cm 2 g" 1 [396J. 

Although different mass estimates agree well in some clusters (e.g. ll86l 12201 123 11 12281 
130711198115181 ). significant discrepancies between cluster masses derived from weak lensing 
and X-ray observations are often found H2741 [6l 12091 13911 11491 and interpreted as signalling 
dynamical processes in unrelaxed cluster cores. Of the 38 clusters in the X-ray selected 
sample studied by 11961 . - 30% show signs of dynamical activity, and more than 50% are 
strong lenses. Based on a sample of 24 clusters between redshifts 0.05 and 0.31, ll95l 
claim that clusters with temperatures < 8 keV show good agreement between different mass 
estimates, while hotter clusters do not. Spectacular examples of massive, distant clusters are 
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Figure 12. The galaxy cluster A 520 (the cosmic train wreck), in which the galaxies are seen 
displaced from the dark matter (lower left panel) and the X-ray gas (lower right panel). In the 
top panel, the lensing signal (blue) and its contour lines are superposed on the X-ray image 
(red) and the cluster galaxies (orange), (from [295 1) 

presented in H28 11 12271 . Despite their youth, the agreement between X-ray and weak- lensing 
mass estimates hints at established clusters. 

Projected cluster mass profiles obtained from weak lensing are often well fit by the 
isothermal profile [144711 or by both the isothermal and the NFW mass profiles [|83l [851 [6]|. 
As data have improved and samples were enlarged, the NFW profile was most often found 
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Figure 13. The galaxy cluster IE 0657-558 (the bullet cluster) whose galaxies and dark matter 
(left panel) are displaced from the X-ray gas (right panel), (from [80|) 



to describe the lensing data very well f[85l 130011 . When based on weak lensing alone, 
NFW concentration parameters tend to be somewhat lower than theoretically expected 
ll86l 12001 12281 130111 . which may be due to intrinsically triaxial cluster halos lf8~6ll . However, 
there is an increasing number of clusters for which NFW profiles with reasonably high 
concentration parameters are deduced (e.g. Il85l [H). It83l find the more massive of six 
high-redshift clusters less concentrated than the less massive ones, which is also expected 
from theory. ll9~7Tl fit the generalised NFW profile to six massive clusters at z — 0.3, 
finding a central double-logarithmic slope a = -0.9 ... - 1.6 at 68% confidence. Assuming 
a = -1, the concentration parameters are well in the expected range, i.e. 5 ... 10 depending on 
cluster mass. The overconcentration problem mentioned earlier persists in many cases when 
constraints from strong and weak gravitational lensing are combined. 

Large-scale structure in front of and behind galaxy clusters is projected onto them and 
can affect weak-lensing mass determinations. Using large-scale structure simulations, [133811 
estimate that weak-lensing mass estimates exceed real cluster masses by several tens of 
per-cents due to the added large-scale structure. H1961 estimated that projected large-scale 
structure approximately doubles the error budget for weak-lensing cluster mass estimates. 
However, cluster mass profiles are affected by cluster substructures and asymmetries only at 
the per-cent level I125T1I81I 

We have seen in the discussion of strong cluster lensing that clusters at moderate and 
high redshifts, z > 0.8, are already remarkably efficient strong lenses. The first weak- 
lensing mass map of a cluster at such high redshift (MS 1054-03 at z = 0.83) was produced 
by H283L The weak-lensing signal of many similarly distant clusters was measured since, 
typically confirming the presence of well-developed, massive and compact clusters at that 
epoch Il84l 11481 12811 13071 12271 . but also frequently indicating violent dynamical activity in 
cluster cores [DM HH |228l . 

Occasionally, detections of clusters with very high mass-to-light ratios (e.g. HI 19111291 1 
are claimed and raise the question whether cluster-sized dark-matter halos may exist which 
are so inefficient in producing stellar or X-ray emission that they are invisible for anything but 
gravitational lensing. [11 1511 detected a peak in the weak-lensing signal 7 arc minutes south 
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of the cluster A 1942 where no optical or infrared emission could be found H157II . A more 
recent analysis 1149 1H revealed a discrepancy between ground- and space-based data which 
remains unresolved. A similarly dark weak-lensing signal peak was discovered next to the 
high-redshift cluster CI 1604+4304 by H482L However, another tangential shear alignment 
potentially revealing a dark halo H344H was meanwhile found to be spurious 1111311 , and a new 
analysis of the A 901/902 supercluster by HI 841 did not confirm the anomalous mass-to-light 
ratio in one of the three matter concentrations found by H158L 

|[422|| introduced the aperture-mass statistic specifically for detecting dark-matter halos 
through their weak-lensing signal. The aperture mass is a weighted integral within a 
circular aperture over the shear component tangentially oriented with respect to the aperture 
centre. When applied to numerical simulations, the aperture-mass statistic turned out to be 
highly efficient in finding group- and cluster-sized halos, although the completeness of the 
resulting halo catalogues has to be balanced against the frequency of spurious detections 
by carefully choosing the signal-to-noise threshold H3981 |499ll . Spurious detections caused 
by projections of large-scale structures are frequent, but can be substantially suppressed 
by optimal filtering H3171 1375H . H506B report the first detection of a galaxy cluster through 
weak lensing, which was confirmed later through its optical signal. (142111 use the aperture- 
mass technique for confirming the weak-lensing signal of clusters found optically in the ESO 
Imaging Survey. Recent applications of cluster-detection techniques based on weak lensing 
H180[|505[|103[|1431l347[|420l have returned catalogues with varying degrees of contamination 
by spurious detections. Their statistical analysis is likely to provide important information on 
the evolution of non-linear cosmic structures in the near future H3161I1041 . 

[PT8l showed that the detection efficiency of the aperture-mass technique varies strongly 
with the density profile of the dark-matter halos, allowing a statistical discrimination between 
isothermal and NFW profiles. [134811 found 4.9 + 2.3 dark-matter halos in a field of 2.1 square 
degrees taken with the Subaru telescope, which is consistent with expectations based on CDM 
models and NFW density profiles [266J . 

As mentioned before, cosmological models with dynamical dark energy cause dark- 
matter halos to be more concentrated compared to models with cosmological constant. While 
this should in principle lead to a higher number of weak-lensing halo detections in dark-energy 
cosmologies and thereby provide a way for discriminating cosmological-constant from dark- 
energy models, the expected sensitivity is very weak due to competing effects j|22ll498H . 1149711 
argue that clusters in formation, which are not virialised yet and thus under-luminous, may 
be detected through weak lensing. They suggest this as an explanation for the potential dark 
clusters found by H115II and [482] and argue that cosmological constraints could be placed by 
comparing the numbers of visible and dark clusters. 

6. Large-scale structures 

Weak gravitational lensing by large-scale structures is covered by several dedicated reviews, 
highlighting different aspects of this rich and quickly developing subject H3201 1271 12051 13221 
13991 12011 I355II . We can only summarise the most important aspects here and refer the 
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interested reader to those reviews for further detail. 
6.1. Expectations and measurements 

Being inhomogeneously distributed in the Universe, matter on scales even larger than galaxy 
clusters must also gravitationally lens background sources. Early studies [1451 13391 12351 
calculated the ellipticities and ellipticity correlations expected to be imprinted on the images 
of background galaxies, and found them to be of order a few per cent on arc-minute angular 
scales. In a first attempt at measuring this weak cosmological lensing signal, 113541 could 
place an upper limit in agreement with theoretical expectations. 

Since weak cosmological lensing is highly sensitive to the non-linear evolution of 
the large-scale structures 1122411 . numerical simulations had to be carried out for precisely 
estimating the amplitude of the signal and the shape of the ellipticity correlation function 
(e.g. B3~l 12251 1 1681 148311 18810 . The cosmological potential of large weak-lensing surveys was 
quickly pointed out ll43l 12361 1485H . emphasising the possibility of measuring in particular the 
matter density parameter Q an d the amplitude <x 8 of the dark-matter power spectrum. 

H429II announced the detection of a coherent shear signal in the field of the radio galaxy 
PKS 1508-05 which they interpreted as being caused by large-scale structure lensing. The 
breakthrough came soon thereafter, when several different groups almost simultaneously 
reported the measurement of the cosmic-shear correlation function IfTTl 14871 15071 [306 ] . 
Given the difficulty of the measurement and the different telescopes, cameras, and analysis 
techniques used, the agreement between these results and their compatibility with theoretical 
expectations was exciting and encouraging. 

Cosmological parameters were soon derived from these first cosmic-shear measurements 
H489B , finding <x 8 > 0.7 and Q ^ 0.4 for spatially-flat cosmological models. Two-point 
statistics of the cosmic shear are approximately proportional to the product ct 8 £2q, i.e. they are 
degenerate in these two parameters. This degeneracy can be lifted using third-order statistics 
such as the skewness H4861 12471 . which arises because the non-linear evolution of cosmic 
structures leads to non-Gaussianity in the weak-lensing signal. Non-Gaussianities were first 
detected by [|4"2l in the Virmos-Descart weak-lensing survey. 

Much effort was devoted to calibrating weak-lensing measurements, to designing 
optimal cosmic-shear estimators and studying their noise properties. HI 1411 used numerical 
simulations to show that relative accuracies of 10% . . . 15% can be reached by cosmic-shear 
measurements. A method for estimating the weak-lensing power spectrum inspired by the 
CMB data analysis was proposed by H214|| . [190114421 1 172H investigated how non-Gaussianity 
can affect cosmological parameter estimates from the cosmic-shear power spectrum. Different 
estimators for the two-point statistics of cosmic shear and their correlation matrices were 
dicussed by Il427l 12341 

Numerous weak-lensing surveys have meanwhile been conducted. A non-exhaustive 
selection of the results obtained on cr 8 for fixed CIq = 0.3 is given in Tab. [2} Although most 
values of cr 8 agree within the error bars, the scatter is still substantial. This is at least partially 
due to remaining systematics in the data analysis, as will now be discussed. 
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Table 2. Non-exhaustive selection of results for eg extracted from weak-lensing surveys. 
Q.q - 0.3 is adopted throughout. erg scales with £2o approximately like erg oc Q Q 5 . 
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Figure 14. Constraints on two cosmological parameters, the normalisation parameter erg and 
the matter-density parameter Qmo, obtained from cosmic-shear correlation functions in two 
different surveys of gravitational lensing (left panel: from [37 1, right panel: from [137]) 



6.2. Systematics 

Numerical simulations show that the Born approximation is valid to very good approximation 
H106L 14441 1188H . While this is so, the effects of weak lensing can be summarised by a scalar 
potential. Then, only such distortion patterns can be caused by weak lensing which can be 
described by derivatives of a scalar potential. This leads to the E and 5-mode decomposition 
discussed in the introduction. Significant B modes in the data are interpreted as remainders of 
undetected or incompletely removed systematics. 

More or less significant B modes have been found in almost all weak-lensing surveys. 
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Until recently, their origin was unclear. [428] showed that source clustering could cause a B- 
mode contribution to the distortion, but not of sufficient strength on arc-minute scales to fully 
account for the observations. H197M showed that incomplete correction for the anisotropics 
in the point-spread function of the imaging system could give rise to a substantial B mode. 
Application of an improved model for the point-spread function to the Virmos-Descart weak- 
lensing survey caused the B mode to disappear [488]. It thus seems that much of the 5-mode 
problem, which was discussed at length in the literature, was due to insufficient correction for 
the distortions imprinted by the imaging system. However, incomplete knowledge of cosmic- 
shear correlation functions, which are necessarily always limited to finite fields, can also lead 
to production of B modes [248 ] . Methods for removing this 5-mode contamination were 
developed by IB351H55I . 

There are (at least) five important sources of systematic error in weak-lensing 
measurements: distortions by the telescope, miscalibrated distortion measurements, 
insufficient knowledge of the non-linear power spectrum, insufficient information on the 
redshift distribution of the background galaxies, and intrinsic galaxy alignments. Extensive 
literature exists on all of these effects. Distortions by the telescope, summarised as 
anisotropics of the point-spread function, are typically corrected by measuring the shapes 
of stars in a data field and fitting the measured distortions by low-order polynomials. More 
recently, H222I showed how the correction for anisotropic point-spread functions can be 
improved. The calibration of shear measurements has been addressed by several large-scale 
efforts, i.e. the STEP programme lfT86l IBTTTl and the GREAT08 challenge (SO. Spatially 
varying calibration errors were addressed by H162II . and [490] studied that sampling errors in 
the redshift distribution have a similar effect as shear calibration errors. Precise knowledge of 
the non-linear power spectrum is important on angular scales smaller than ~ 10' H188L where 
baryons can also change the dark-matter power spectrum noticeably (123311 . 

Analyses of cosmic-shear measurements assume that intrinsic galaxy ellipticities are 
uncorrelated such that they average to zero when several images are combined. However, 
galaxies being physically close to each other are also expected to have their shapes aligned, 
e.g. by the tidal field of the large-scale matter distribution into which they are embedded. The 
potential effect of intrinsic rather than lensing-induced galaxy alignments depends obviously 
on the depth of the survey. Deep surveys project galaxy images along light paths which are 
substantially longer than any large-scale structure correlation scale and thus suppress any 
spurious signal due to intrinsinc alignments of physically neighbouring galaxies. In shallow 
surveys, however, intrinsic source alignments may substantially contaminate any weak-shear 
signal ||9~3l 11761 1701 |92l . Measurements of intrinsic alignments were found to agree well 
with these predictions H2701 1621 . Recent numerical H232B and analytic studies H269H claim 
strong intrinsic alignments of galaxy halos. Comparing simulations and data, [112211 specified 
that theoretical predictions based on the orientations of dark-matter haloes overpredict the 
alignment by a factor of ~ 2, but agree well with measurements if the central orientations of 
galaxies placed in these haloes are being used instead. 

Possibilities for removing the signal contamination due to intrinsic alignments were 
discussed extensively. They advocate using photometric redshifts to remove physically close 
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pairs of source galaxies from the analysis 112491 11851 . An application of this technique to 
the multi-colour Combo- 17 survey [ 1821 showed intrinsic alignments near the lower end of 
the theoretical predictions and slightly lowered the value of <x 8 previously derived from these 
data [[611 to 0.67 + 0.1 for Q, = 0.3. [ 1 95 1 point out that foreground galaxies are aligned 
with the large-scale structures lensing background galaxies, thus giving rise to a higher- 
order alignment between galaxies at different redshifts. They quantified this systematic effect 
H2981I191I and showed that it can substantially reduce the measured shear signal, leading to a 
likely underestimate of the <x 8 parameter by several per cent. 

6.3. Perspectives 

For several years now, cosmological parameters are obtained with high accuracy mainly 
from combined analyses of CMB data, large-scale galaxy surveys and the dimming of 
type-la supernovae with redshift. What is the role of weak gravitational lensing in this 
context? It should be emphasised that parameter constraints from the CMB alone suffer from 
degeneracies which can only be broken using additional information. Measuring directly 
the dark- matter density and the normalisation of its fluctuations, gravitational lensing adds 
constraints which substantially narrow the parameter ranges allowed by the CMB, as H213I 
pointed out. [[891 combined weak-lensing data from the Red-sequence Cluster Survey with 
CMB data from the WMAP satellite to find cr 8 = 0.89 + 0.05, Q = 0.30 + 0.03 and a Hubble 
constant of Ho = (70 + 3)kms _1 Mpc" 1 . H473H predicted the accuracy of joint cosmological 
parameter estimates using CMB data in combination with the weak-lensing constraints from 
the CFHT Legacy Survey. 

The exploitation of higher than second-order statistics will become increasingly 
important for breaking degeneracies in the weak-lensing parameter estimates and the 
investigation of non-Gaussianity developing due to non-linear structure growth H517[ 14711 
I108H . Much valuable information must also be contained in the morphology of the two- 
dimensional weak-lensing pattern. 

Among the most exciting perspectives of weak lensing is its potential to study the three- 
dimensional distribution of dark structures. Although lensing observables measure the two- 
dimensional, projected density distribution, selecting sources at different distances allows 
structures along the line-of-sight to be resolved. Sufficiently accurate information on the 
distance to the source galaxies is provided by photometric redshifts. H454II show that the 
accuracy of cosmic-shear parameter estimates can be improved by a factor of 5 ... 10 by 
splitting the source galaxies into only four redshift bins. H175il developed an elegant formalism 
for extracting three-dimensional information from weak-lensing data (see also H38510 . and 
H472II recovered the three-dimensional matter distribution in the field of the Combo- 17 survey, 
using photometric redshifts as distance indicators for the source galaxies. 

The possibility to extract three-dimensional information from weak-lensing data opens a 
way for studying the growth of cosmic structures along the line-of-sight from the distant and 
past universe. This, in turn, is sensitive to the change of the dark-energy density with time. 
Thus, sufficiently accurate and wide-field weak-lensing surveys will allow the dark energy 
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to be constrained H215L The sensitivity of such experiments is very promising 0561 [341 . in 
particular if higher-order statistics are considered. This is one of the strongest motivations for 
a weak- lensing survey from space H314[ 14011 . The impact of several sources of systematic 
error on the determination of dark-energy properties is of great interest in this context 
H2501 11771 1571 125311 . Enormous progress is expected regarding the control of systematics 
and the suppression of noise from planned space-based surveys of large fractions of the sky, 
such as the survey to be undertaken by the proposed DUNE-Euclid satellite mission [140011 . 

Planned wide-area surveys of the sky in radio wavebands are expected to yield 
information on both cosmology and the distribution of neutral hydrogen during reionisation 
[386|. Weak lensing in the data of the FIRST radio survey undertaken at 20 cm wavelength 
was recently detected by IT76ll . 

6.4. Cosmic magnification 

Besides distorting the images of distant galaxies, large-scale structures also magnify 
background sources. To linear order, the power spectrum of cosmic magnification is simply 
four times the power spectrum of the cosmic shear, i.e. both contain the same amount of 
cosmological information. Cosmic shear, however, is much more easily measurable than 
cosmic magnification because it can reasonably be assumed that the ellipticities of distant 
galaxies average to zero. The intrinsic flux of any given source being unknown, cosmic 
magnification is much harder to identify. 

Currently the most promising (and perhaps the only) method for detecting cosmic 
magnification is the magnification bias. If a population of distant sources, e.g. quasars, is 
observed within a solid angle 6Q, on the sky where the magnification is /u, fainter sources 
become visible there, but their number density is reduced because the solid angle is stretched 
by the magnification. The net effect depends on how many more sources the magnification 
lifts above the flux threshold of the observation. If the number-count function of the sources 
is sufficiently steep, the dilution is outweighed and the magnification causes more sources to 
be visible, while sources with flat number counts are diluted and appear less numerous than 
without the magnification. 

The spatial galaxy distribution follows the dark-matter structures which act as lenses on 
background sources. Bright quasars, for example, have a steep number-count function and 
thus appear more numerous behind magnifying dark-matter overdensities. The correlation 
of foreground galaxies with the lensing matter then leads to a positive cross -correlation on 
angular scales of arc minutes and larger between distant QSOs and foreground galaxies which 
are physically separated by cosmological distances. 

The cross-correlation function between background sources and foreground galaxies was 
derived in linear approximation by [14J and H 10711 : see also H161H . However, magnification 
is non-linear in shear and convergence. Computing the theoretical expectation of the 
cosmic magnification accurately is thus considerably more complicated than for the cosmic 
shear. [287J derived second-order corrections to the linear cosmic-magnification statistics 
and found that the linear approximation was systematically low by 20% . . . 30%, which 



Gravitational Lensing 



56 



was confirmed by H470B in a fully non-linear treatment based on the halo model of the 
dark-matter distribution. The accuracy of analytic calculations can be calibrated using 
numerical simulations of light propagation through large-scale structures (e.g. [fT2l 128711 ). 
Sufficiently accurate analytic and numerical calculations of cosmic magnification yield typical 
magnifications of < 10% for point sources at redhifts z s — 1. While most emphasis was 
put on cross-correlations between QSOs and galaxies, cosmic magnification may also induce 
detectable cross-correlations between background and foreground galaxies [13491 13501 . H2211 
discussed survey strategies for detecting cosmic magnification. 

11241 confirmed that the QSO-galaxy correlations detected by H437H and H138B can indeed 
be explained in terms of gravitational lensing by large-scale structures. Subsequent more 
detailed correlation studies confirmed that the correlations showed the significance variations 
expected from the lensing hypothesis 11251 , and revealed cross-correlations of distant QSOs 
with infrared [1261 1321 and X-ray emission Il28l . While the existence of QSO-galaxy cross- 
correlations was thus firmly established, their amplitude and angular scale was typically 
found to be substantially too high [50Ql 13671 l36l 11451 . although some analyses obtained 
the theoretically expected results H408[ 13661 . The observational situation was thus utterly 
confused. HI 341 found evidence for cosmic shear in the vicinity of five distant QSOs. Dust 
absorption in foreground galaxies may be responsible for similar anti-correlations as expected 
from weak lensing of faint background QSOs (e.g. [|94l ), and selection effects may cause 
correlations as well as anti-correlations (e.g. H128I0 . The detailed relation between dark- 
matter halos and their galaxy content adds further uncertainty H223H . 

Due to the weakness of the signal and the possibility of confusing it with other 
effects, large multi-colour surveys such as the SDSS were expected to be necessary for an 
unambiguous detection of cosmic magnification and its exploitation in terms of cosmological 
parameters OH I285H . The degeneracy between the density and bias parameters Q, and b 
can be broken determining the three-point cross-correlation between QSOs and galaxy pairs 
[286J. Clear and unambiguous evidence for cosmic magnification was finally detected by 
[143 111 who measured the cross-correlation between a colour-selected sample of - 200, 000 
distant QSOs and foreground galaxies in 3, 800 square degrees of the SDSS data and found 
a signal significant at the 8-<x level and in complete agreement with theoretical expectations. 

A very promising manifestation of cosmic magnification was recently discussed by 
[288J. Distant QSOs magnified by intervening matter are more likely to show absorption 
features in their spectra due to the gas associated with the lensing structures. Using the 
QSO sample from the 2dF survey, H288H demonstrated evidence for lensing by Mgll and 
Fell absorbers along the lines-of-sight to the QSOs. This was further studied by B23B and 
discovered in the SDSS data by 0241 . 

H4841 pointed out that cosmic shear and magnification provide complementary 
cosmological information. Using deep data of the CFHLS, [189] detected the effect of cosmic 
magnification by cross-correlating distant, so-called Lyman-break galaxies with foreground 
galaxies. 
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6.5. Gravitational lensing of the CMB 

This subject has been comprehensively reviewed by 1127311 . Pioneering studies 11721 I7T1 
14381 14401 13101 pointed out that the CMB is expected to be weakly gravitationally lensed 
at a measureable level and developed methods for computing the effect of lensing on the 
CMB power spectra. Lensing of the CMB can be seen as a diffusion process whereby 
structures are blurred, but also created on small angular scales in the Silk-damping tail 
H336L The effect of lensing on the CMB cannot be identified in the CMB power spectrum, 
but in higher-order statistical measures 11391 |40l . Distortions of the probability distribution 
function of CMB temperature fluctuations, higher-order correlations, quadratic filtering and 
maximum-likelihood approaches were developed to identify lensing in the CMB at the 
statistical level and in two-dimensional reconstructions [11641 12121 12461 1194[ I193H . It was 
pointed out that lensing of the CMB must be treated in spherical geometry on the full sky 
1121 1 L 1741 1373[ 1731 if percent-level accuracy is to be achieved. Since lensing modifies the 
CMB power spectrum from which cosmological parameters are derived at high precision, 
undiscovered or uncorrected lensing leads to biases in cosmological parameters H272L On 
the other hand, lensing can break cosmological parameter degeneracies in the CMB data 
H465L Gravitational lensing of the CMB by fully non-linear cosmological structures was 
simulated on the full sky by [|68l |67ll , who showed how even large-scale 5-mode CMB 
polarisation can be created by small-scale non-linear structures. Simulations were also used 
to study how well CMB lensing can be recovered in presence of systematic effects [138711 . So 
far, there are no published direct detections of CMB lensing, but there is indirect evidence 
with significance around 3-cr obtained by cross-correlating the CMB with different samples 
of distant foreground sources used for tracing the potentially lensing large-scale structure 
[460, 190]. Direct evidence for CMB lensing should soon be discovered in the data of the 
Planck CMB mission. 

7. Summary 

Many are the applications of gravitational lensing to cosmology, and the results are numerous, 
as the preceding discussion has shown. A review like this must be based on a subjective 
selection which is necessarily biased to some degree. Within these limitations, I summarise 
the results as follows: 

• Microlensing experiments in the Galaxy have shown that, although massive compact 
objects exist in its halo, they are insufficient to make up all the dark matter in the Galactic 
halo. These studies have extended towards the Andromeda galaxy M 31. It is not clear 
yet what fraction of the observed microlensing must be attributed to self-lensing by the 
visible stars. Low-mass planets have been detected by means of microlensing. 

• Central density profiles of lensing galaxies are well described as isothermal within the 
radial range where they produce multiple images. Their cores are thus more concentrated 
than CDM predicts. This indicates that galaxy density profiles have been steepened by 
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baryonic physics. At larger radii, weak galaxy-galaxy lensing shows that the isothermal 
density profiles steepen and approach the NFW density profile shape. 

• Galaxy-galaxy lensing finds large halo sizes with radii of > 200 h~ l kpc. Halos of cluster 
galaxies seem to be smaller, as expected. The biasing of galaxies relative to the dark- 
matter distribution is found by galaxy-galaxy lensing to be almost scale-independent, or 
gently increasing with scale. 

• Galaxies have to be structured in order to explain multiple-image geometries and the high 
fraction of quadruple compared to double images. Anomalous flux ratios of quadruple 
images seem to be best explained by lensing, but simulations show that the expected level 
of substructure is insufficient to explain the observed anomalies. 

• Measured time delays between multiple images lead to an interesting conflict between 
the lensing mass distribution and the Hubble constant: Isothermal profiles yield Hubble 
constants which are substantially too low, and lens models giving compatible Hubble 
constants have too steep mass profiles. It seems that this conflict can be resolved allowing 
perturbations of the density profiles. 

• The statistics of distant sources multiply imaged by galaxies is sensitive to the 
cosmological parameters. Recent applications of this method showed agreement with 
a low-density universe with cosmological constant. 

• Galaxy clusters have to be asymmetric, and they must be dominated by dark matter which 
is more broadly distributed than the cluster light. Cores in the dark-matter distribution 
must be small or absent. Frequent and substantial discrepancies between lensing and X- 
ray mass determinations are most likely signalling violent dynamical activity in clusters. 

• It seems that galaxy clusters in the "concordance", low-density spatially-flat 
cosmological models cannot explain the observed abundance of gravitational arcs. 
Clusters need to be highly substructured and asymmetric, and their dynamics temporarily 
boosts their strong-lensing efficiency. Yet, theoretical expectations fall substantially 
below extrapolations of the observed number of arcs. Surprisingly massive and compact 
clusters which are significant weak and powerful strong lenses exist at redshifts z - 0.8 
and above. 

• Although cluster density profiles inferred from strong and weak lensing do typically not 
contradict expectations from CDM, isothermal density profiles are not ruled out by strong 
gravitational lensing. Claims of flat central profiles are not supported by reasonably 
asymmetric models. 

• Typical mass-to-light ratios derived from weak cluster lensing range around 200 in 
solar units, but very high values have occasionally been found. While this may indicate 
a separation of gas from dark matter in cluster mergers, the possible existence of dark 
clusters is intriguing. 

• Cosmic shear, i.e. the distortion of background-galaxy images due to weak gravitational 
lensing by large-scale structures, has been detected and found to be in remarkable 
agreement with theoretical expectations. It has enabled constraints on the matter-density 
parameter and the normalisation parameter <x 8 of the dark-matter fluctuations. 
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• Systematic effects such as image distortions in the telescope, calibration errors on shape 
measurements, insufficient knowledge of the non-linear matter power spectrum and the 
redshift distribution of background galaxies and possible intrinsic alignments of source 
galaxies are important and substantial and need to be carefully corrected. 

• Joint analyses of weak lensing and CMB data allow parameter degeneracies in both types 
of experiment to be lifted. When combined with photometric redshifts of source galaxies, 
three-dimensional reconstructions of the large-scale matter distribution become possible. 
This will also allow constraints on the dark energy. 

• Cosmic magnification, which is more complicated to measure than cosmic shear, can be 
quantified by the magnification bias. It has been detected, and most recent measurements 
are also in excellent agreement with theoretical expectations. 

• Gravitational lensing of the Cosmic Microwave Background is inevitable and affects 
cosmological parameter estimates obtained from the CMB at the per cent level, if 
uncorrected. It broadens the peaks in the CMB power spectra, creates small-scale 
temperature fluctuations in the Silk damping tail and converts part of the it-mode 
polarisation of the CMB into B modes. CMB lensing has been marginally detected at 
the expected level by cross-correlating the CMB with distant foreground sources. 
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Appendix A. Light deflection in Newtonian approximation 

Usually, the equations of gravitational lensing are derived from the deflection angle of a point 
mass in Newtonian approximation. To establish the link between this derivation and that based 
on the equation of geodesic deviation given above, we show here two alternative derivations 
for the Newtonian deflection angle. One is based on the geodesic equation, the other on 
Fermat's principle. 



Appendix A. 1. Geodesic equation 

The geodesic equation asserts that the null wave vector k is parallel transported along the light 
ray, 

V k k = => k(k a ) + kPafpik) = . (A.l) 



For evaluating it in the Newtonian metric (28), we conveniently return from the dual basis 
(29) to the coordinate basis. Then, k = oj(1, e), and we find 

(d^ + e'd^j co + coe j (pj = , 

(d„ + e l d) e j = - 2 - eXe'fr)] (A.2) 
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for the time and space components of the geodesic equation, respectively. The second 



Eq. (A. 2) can be re- written in the conventional way 

d/ = -2ff<f> . (A3) 

Identifying the tangent vector e l with 
dx' 

e> = (A.4) 
drj 



Eq. (A3) already reproduces Eq. (40). 



Appendix A.2. Fermat's principle 

A static metric can be written in the form 

ds 2 = #ooC 2 d? 2 + g ij dx'dx i , (A.5) 

where the g Q0 and g t j are independent of time. Let x M (A) be a light path parameterised by the 
affine parameter A and i p = dx^/dA, then 

(x, x) = (A.6) 

and the light paths are given by the variational principle 

"viz 

5 1 (x,x)dA = 0. (A.7) 



(x, x) 



In this stationary space-time, we can choose the affine parameter such that 

9^ = -1 . (A.8) 
dA 

If we carry out the variational principle giving up the condition that 6t = at the end points 
of the paths, 



r A 2 i 2 rh. 

= 6 (x,x)dA = -26t = -26 dt . (A.9) 

Since 

- g QC 2 dt 2 = gijdx'dx j = dcr 2 (A. 10) 

for light rays, Eq. ( |A.5| ) implies 



^ do- 
V _ #oo 

This gives Fermat's principle for a static space-time in general relativity. For the Newtonian 



6 f -^= = . (A.ll) 



metric Eq. (28) 



-goo =1+20, dcr = (l -2(f>) l/2 \iv\ dA , (A.12) 

such that 

5 f -\r—^r\^\ dA = s f n(w)\iv\dA = , (A.13) 
J M y 1 + 2<p J M 
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where 

m =i^^- 2 * <a - i4) 

is the effective index of refraction. The Euler-Lagrange equation of the variational principle 



n— = Vn - (e ■ Vn)e = V ± n (A.15) 



( |AJ3] >is 

dA 

if we scale the affine parameter A such that w = e has unit length, e 1 = 1 . Writing 

iv ± n = V ± lnn = V ± ln(l -20) * -2V ± , (A.16) 

we reproduce the equation of motion 

dJ=-2V ± <p (A.17) 

derived above directly from the geodesic equation. Note that, to first order in 0, we can 



identify de/dA with de/drj; see Eq. (A. 3 ). 



Appendix A3, Deflection angle of a point mass 
For a point mass M at the coordinate origin, 

c z r 2r 

with the Schwarzschild radius R s . Anticipating a small deflection, we can use Born's 
approximation, integrate V ± cp along the unperturbed, straight light path and thus write 



Eqs. ( |A.3] ) or flA.17| ) in the form 

ff. / r \ 

(A. 19) 



de Rj x\ R 



dz r 3 \y ) ( x 2 + ^2 + z 2)3/2 ^ y 

if we introduce the spatial coordinate system such that the light ray propagates into the e z 
direction. Integration over z gives the deflection angle 

d = -^b, b = ( X ). (A.20) 



b 2 ' \y 

Since the field equation is linear in the Newtonian limit, the deflection angle of more 
complicated mass distributions is a linear superposition of point-mass deflection angles. 
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